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Abstract 

I review a Constituent-Quark-Meson model (CQM) for heavy meson decays, outlin- 
ing its characteristics and the calculation techniques developed for it. The strength 
of this effective model, is that it enables to evaluate heavy meson decay amplitudes 
through diagrams where the heavy mesons are attached at the ends of loops con- 
taining heavy and light quark internal lines. The phenomenological applications are 
presented in detail, trying to give a self-contained operative picture of the model. 
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1 Introduction 



During recent years, heavy meson physics has received a wide attention both from theory and 
experiment. This is because it helps the comprehension of many open problems of the standard 
model and can also act as a passage in the domain of new physics. Many experiments on B 
physics already at work or near to be started, BaBar, Belle, CLEO III, Hera-B, CDF-DO and 
those planned to begin after 2005, ATLAS, CMS, LHCB and BTeV confirm this interest 0]. 
B physics has had an important role also in LEP I that has registered about 10^ — > bb 
events B decays offer the framework for investigating in detail the field of CP violations 
and for determining CKM (Cabibbo-Kobayashi-Maskawa) matrix elements. In particular, rare 
B decays, those in which there is no charm in the final state, are relevant for the research 
of signals of new physics In fact, the Standard Model predicts that rare B decays (the 
Cabibbo suppressed or the penguin induced decays) should be strongly suppressed, therefore, 
any anomalous increasing of branching ratios could be due, for example, to the existence of 
new particles, external to the standard model spectrum because interacting at higher energy 
scales. 

The amplitudes governing heavy meson decays are theoretically calculated mainly using 
lattice QCD methods and the SVZ (Shifman-Vainshtein-Zakharov) sum rules Q . 

The lattice QCD program [||, is that of computing the QCD partition functional summing 
over a representative ensamble of gauge fields and fermionic field configurations; the action is 
written in discrete form modelling the entire space-time as a four-dimensional grid where the 
distance between nearest neighboring sites is a and the linear dimension is L ~ il4 , 17 being the 
four volume of the grid. In principle, considering a sufficiently large number of configurations 
and simulating a very close (a — > 0) and large {L oo) grid on a calculator, amounts to build 
a calculation framework nearly resembling that of continuous QCD. In practice, there are many 
technical problems: some of them have to do with computer power, some with the continuous 
limit of the results obtained on a discrete space-time grid. 

In the ordinary hadronic matter, the quarks are not very far from each other, therefore, in 
ordinary circumstances, it is not essential to consider the complex QCD dynamics giving rise 
to the Abrikosov chromoelectric flux tubes thought to be responsible for quark confinement. In 
this situation valence quarks are weakly interacting with QCD vacuum fluctuations. The SVZ 
method aims at determining the parameters and the regularity of ordinary mesons and baryons 
through an expansion of the correlation functions, written in terms of dispersion intergrals, in 
a power series controlled by the parameter (the strong coupling constant), plus power cor- 
rections expressed through the vacuum condensates {Cj^^, qq, qaGq, ..). It is believed that the 
vacuum condensates contain the most relevant non perturbative effects of the QCD vacuum. 
Invoking the concept of parton-hadron duality, this expansion must be compared to the phe- 
nomenological expressions for the correlation functions. It is this comparison that allows to 
extract quantitative information on 2, 3, ..-points correlators, i.e., on all possible observables. 
One of the main drawbacks of SVZ sum rules is the difficulty one meets in computing the 
theoretical error due to the ambiguous choice concerning the truncation point of the series 
expansion. 

This work is devoted to introduce an effective Constituent-Quark-Meson model based on a 
Lagrangian incorporating the symmetries of heavy quark effective theory, the chiral symmetry 
in the light quark sector, see section ^, and, as is discussed in section ^, dynamical information 
derived from an underlying Nambu-Jona-Lasinio interaction. In section ^, together with the 
discussion of calculation techniques used for computing some relevant loop-integrals, it is shown 
how the determination of strong coupling constants, parameterizing the low energy effective 
hadron Lagrangian, proceeds through a comparison of the low energy matrix elements with the 



2 



CQM computed amplitudes: CQM plays the role of a fundamental model (since it contains, 
besides meson fields, also the elementary heavy and light quark fields) with which the hadron 
theory must match at higher energy, see discussion in section |2.1.1j . With respect to lattice 
QCD and SVZ sum rules, CQM is a rough approach that, anyway, has shown to be a quite 
reliable and easy-to-use method. 

One of the very common problems of quark models Q, is that of associating theoretical 
errors to predictions. This topic is discussed for CQM in section ^, together with the problem 
of defining the light constituent quark mass. The constituent quark mass is typically heavier 
than the current mass, appearing in the QCD Lagrangian (and related to the Higgs field VEV): 
one can think of a constituent quark as of a current (bare) quark dressed by a cloud of virtual 
particles generated by strong interactions Q. The mechanism dressing the bare quark and 
giving the constituent quark its mass value, is an intrinsic feature of the model itself. 

Section ^ is devoted to the study of exclusive semileptonic decays of B mesons through 
the CQM model. Here are examined processes involving b ciu and b uiv transitions, 
the former being related to Vcb, the latter to Vub- In particular, CQM has allowed to obtain a 
prediction for the branching ratio of the semileptonic process i? — > ai. 

All existing evaluations of exclusive semileptonic B decays are strongly model-dependent 
or are affected by problems related to the estimation of the theoretical error. Anyway an 
agreement among diverse models, e.g., on the determination of a particular form factor, gives 
rise to a theoretical platform useful for a comparison with experimental data. This could also 
be an alternative approach to the study of rare B decays, considering that the most commonly 
used method to extract Vub through a comparison with data, is the so called end-point-method, 
see, e.g., [^. The idea of the end-point-method is that of eliminating the background due to 
b civ decays while examining the inclusive leptonic spectrum -^^^{b — > uii') in the Eg region 
where the invariant mass Mx of the hadron system emerging from the decay is such to avoid 
decays in a charmed final state: Mx < Md. But, in this region of the energy spectrum, one 
meets technical difficulties related to the Wilson expansion of ^^: one can only compute the 
first terms of this expansion. Higher order terms depend on matrix elements of local operators 
having higher dimensionality, and can at most be estimated by phenomenological models. It 
is possible to show that, in the proximity of the end-point, i.e., in the proximity of a certain 
critical value Mx,c, ah terms in the Wilson expansion are equally important and, for even 
higher values of Ei, the decay cannot anymore be analyzed by Operator-Product-Expansion. 
In the experiments devoted to the determination of Vub^ a kinematic cut on Mx, very near to 
the critical value Mx,c, is used. This means that the determination of Vub is model-dependent 
since it is necessary to be able to estimate the terms having higher dimensionality in the Wilson 
expansion. To avoid this problem, one could think of enlarging the Ei region experimentally 
examined. This could give the possibility of being far from Mx,c, but the price to pay is that 
of a strong growth of the background of events containing charm in the final state. 

CQM gives the possibility of further investigating the exclusive channels — > />, i? — > ai 
and i? — > vr in such a way to enlarge and give more solidity to the platform of model-dependent 
results I mentioned before. 

2 Introduction to the formalism 

2.1 Effective theories 

In this section I will discuss briefiy the general topic of effective theories in particle physics 
with the aim of introducing the basic ideas and tools of the CQM model in the subsequent 
sections. 
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When one calculates the energy levels for an hydrogen atom, the problem to face is that of 
solving the Schrodinger equation for an electron moving in the Coulomb potential generated 
by the positive proton charge: it is not relevant to take account of the inner quark structure of 
the proton. The low energy dynamics of the hydrogen atom does not depend in any relevant 
way on the high energy finer details of the proton inner structure. The proton can be simply 
considered as the static source of Coulombic potential and, in a first approximation, we can 
ignore also its spin and magnetic moment. Doing in such a way, the problem of determining 
hydrogen energy levels presents essentially only one energy scale rrie (the electron mass) and 
the dimensionless fine structure constant a: we have separated out higher energy scales. This 
can be done essentially because of the large separation of the energy scales that usually enter 
into a physical problem. A physical system in which there are different but close to each other 
energy scales, cannot be treated in the same way because even small perturbations can allow 
the system to explore all these scales with similar probabilities. 

A finer calculation of the hydrogen energy levels requires to include in the calculation the 
effect of the spin and of the magnetic moment of the proton. These details are responsible of 
the well known hyperfine structure of the energy levels. We can state that the energy levels of 
the hydrogen atom can be computed ignoring the dynamics acting at scales larger than A, with 
A >> rriea, with an error of order mga/A. The more the desired precision, the higher is A, the 
smaller is the error one makes ignoring the high energy (> A) dynamics. For example, parity 
violation effects at the atomic level are very small since the weak interaction energy scale is 
Mvi/, extremely larger than the atomic energy scale. 

Effective theories ||9|-|l^] are those models conceived to describe the physics of a certain 
system at the energy scale of the experiment through which one studies it, i.e., at the level of 
accuracy chosen to experimentally examine the system. In this sense, the atomic physics of 
the hydrogen atom is an effective theory of the hydrogen. 

Effective models succeed in giving reliable phenomenological predictions where fundamental 
theories have many more technical and sometimes principle problems. Quantum-Chromo- 
Dynamics (QCD) is the most important example of a fundamental theory, i.e., a theory derived 
from first principles, describing the intimate nature of strong interactions and the building fields 
of matter, that has deep troubles in dialing with the low energy hadron world. This is due 
to the still partial theoretical comprehension of the confinement mechanism of quarks in the 
hadronic matter. Therefore, to deal with hadrons, it is necessary to implement some low energy 
model, effective in the energy regions where the hadronic processes one wants to study take 
place. 

A low energy effective theory of hadrons is anyway a relative of QCD, since it incorporates 
the symmetry properties required by the fundamental theory. The hadron effective Lagrangian 
must therefore be Lorentz invariant, the S-matrix must be unitary, the VCT symmetry must be 
obeyed and it has to show chiral symmetry in the limit in which light current masses are sent 
to zero. New symmetry properties could also emerge in the effective theory being absent in the 
fundamental one: the example relevant for this work is that of Heavy-Quark-Effective-Theory 
(HQET), to which is devoted the next section. 

Symmetry properties select an infinite class of Lagrangian interaction terms, only a finite 
number of them being renormalizable. The requisite of renormalizability, crucial for a funda- 
mental theory, is lost in the effective theory approach. 

The origin of non-renormalizable interactions is due to the absence of heavy particles from 
the spectrum of the effective theory. An example comes from Fermi's /3-decay theory, where 
a non-renormalizable four fermion contact term, distorts the high energy interaction mediated 
by the W particle, absent in Fermi's theory. Anyway Fermi's theory works extremely well 
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at the energy scales of nuclear processes. The masses M of the heavy particles, excluded 
by the effective theory spectrum, may appear as energy cutoffs A = M suppressing the non 
renormalizable terms by factors of E/M, E being the characteristic low energy scale of the 
processes described by the effective theory. For example, the typical energy scale of Quantum- 
Electrodynamics (QED) processes is of order of mg, that is a sufficiently small number to 
explain way QED can be very well considered as a fundamental, renormalizable theory of 
electrodynamic interactions. 

In general terms one can associate to each mass of a known particle a boundary between 
two different effective theories: the anomalous breaking of scale invariance, manifested in the 
peculiar distribution of particle mass values, gives then rise to a tower of separate effective 
theories. For energies below a certain boundary value, one can construct a low energy effective 
theory in which all the particle states above the boundary threshold are excluded from the 
spectrum. Of course, the coupling constants in the interaction terms related to the light fields 
should vary with continuity at the boundaries. 

As we go down in the energy ladder, we meet effective theories containing less fields and 
a larger number of non-renormalizable terms while, in the opposite direction, we find that the 
non-renormalizable terms are progressively more important (less suppressed by E/M) and dis- 
appear at boundaries, where they are substituted by new renormalizable interaction terms. The 
important point is that what happens at high energies doesn't affect the low energy behaviour. 
This picture is deeply explained in |12], [13|. 

The renormalization group method |15] allows to bridge between two effective theories. 
The aim is that of calculating the low energy parameters through the high energy ones. These 
calculations can be explicitly performed only once the high energy theory is weakly coupled. 
The QCD case is therefore complicated because the renormalization group method doesn't allow 
to bridge continuously from the fundamental theory, the QCD, to the hadron effective theories. 
This is why, many times, the hadron low energy effective theory parameters are determined 
by a matching with some other more fundamental model, i.e., some model containing in its 
spectrum also the higher energy elementary particles. These models are not necessarily QCD 
derived, like lattice-QCD or SVZ sum rules. In many cases these models contain hypotheses 
in conflict with the QCD structure. Object of this work is to introduce one of these effective 
models. 

What is important to focus on, is that the proliferation of non-renormalizable terms (the 
irrelevant terms in the modern language), doesn't spoil the predictive power of the effective 
theory. On the contrary, non-renormalizable terms can help in determining the predictive 
power at disposal. 

Here follows an example of how the effective theory approach could make things very easy 
with respect to a fundamental theory approach. 



2.1.1 Photon-photon scattering 

Let us suppose to be interested in understanding how the cross section for the photon-photon 
scattering scales with the energy of the photon in the limit in which this is lower than the 
rest energy of the electron. From an effective field theory point of view, this means that we 
are interested in building an effective theory in which the electron is excluded by the particle 
spectrum. The electron mass acts as the cutoff A = discussed before. 

We therefore only need an interaction Lagrangian containing four photon fields. The sym- 
metry principles instructing us about how to build this low energy effective theory are: Lorentz 
invariance, gauge invariance and the V,C,T symmetries. To the lowest order we can therefore 
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(a) 



(b) 



Figure 1: QED photon-photon (logarithmic divergent) scattering diagram in (a) 
effective interaction in (b). 



The Euler-Heisenberg 



write the so called Euler-Heisenberg Lagrangian: 



mi 



2 r 



ai{F^,Fn^ + a2{F^,Fn^ + 0(l/m: 



(1) 



written in such a way to have the correct mass dimension. ai^2 are the constants multiplying 
respectively the scalar squared and the pseudoscalar squared terms. The presence of F^^F^'^ 
explains why we cannot have odd powers in higher order terms. Every gradient of the pho- 
ton field, in the lowest order Euler-Heisenberg Lagrangian, produces a factor of E^. We can 
therefore argue that the cross section scales as: 



a oc 



1 

7 



(2) 



i.e., a (X The phase space factor 



is needed because a has dimensions of a surface and E, 
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is the only dimensional parameter in the effective theory. The effective approach description is 
the one given in fig. ||(b). Anyway we could calculate the photon-photon scattering in QED, 
see fig. |l|(a), by the virtual electron box (we should add five more diagrams renormalizing its 
logarithmic divergence), i.e., we could calculate the a of the process at high energy using the 
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fundamental theory and then consider the low energy limit of the result. In such a way, through 
the matching of high and low energy Green's functions, we could obtain the constants ai and 
02- If we didn't knew QED, we should have fixed ai and 02 by matching with the experiment. 
In the case of effective hadron Lagrangians, this is the problem: one cannot calculate the 
couplings (like ai and 02), essential for the phenomenological predictions, by a matching with 
QCD. This is why one tries to build effective models, intermediate in energy between QCD 
and the hadron world, that could allow this matching. 

2.2 Heavy quark effective theory 



The physics we are interested in, is that of mesons containing one heavy quark (b or c) |16|- 
p^]. These states present a large separation of energy scales: on one hand we have the heavy 
quark mass niQ and, on the other hand, we have Aqcdj the asymptotic freedom scale, which 
limits the boundary between the strong coupling and the weak coupling region. The heavy 
quark is surrounded by a cloud of light quark states interacting with it through soft gluons 
having momenta of order Aqcd — 200 — 300 MeV. Being mg >> Aqcd, we understand that 
the exchanged soft gluons can resolve only larger distances than the heavy quark Compton 
wavelenght. This means that light quark degrees of freedom are blind to the heavy quark 
flavor (i.e., to mass) and spin. For the light degrees of freedom, the heavy quark is simply a 
static chromoelectric source as the proton is simply a source of Coulombic potential for the 
electron in the hydrogen atom (chromomagnetic effects are suppressed by a factor of l/mg; 
spin-spin coupling terms between light and heavy degrees of freedom are also 1 /mq terms) . 

We can therefore state that light degrees of freedom in an heavy meson have a new symmetry 
property, not remnant of the underlying QCD description, with respect to flavor and spin 
rotations of the heavy quark with which they interact. In particular, the excitation spectra of 
two heavy mesons containing two different heavy quarks Qi and Q2 with mQ^,mQ^ » Aqcd, 
are the same once one overlaps the ground states. Due to flavor symmetry, it happens something 
like the atomic physics independence of the electron structure on the neutron number contained 
in the nucleus. Due to spin symmetry, each excitation level will be a doublet, degenerate in 
the total spin (if light degrees of freedom are not carrying zero spin). 



In a seminal paper by H. Georgi |23], the initial ideas on heavy mesons and flavor-spin 
symmetries |^, are translated in the effective field theory language. The aim is that of 
building a low energy theory where the heavy quark mass is considered infinite, mq 00, 
since it is greater than all other energy scales appearing in the effective theory, while the heavy 
quark velocity, which is practically the same of that of the entire hadron, is a conserved contant 
of motion. Consider: 

p/^ = mnv^, (3) 

the momentum of a meson having mass mn and containing an heavy quark of mass mq. In 
the ruQ 00 limit, mn = mq. Of course, in physical situations, the infinite mass limit is 
not rigorously fulfilled and niH ^ mq . If we suppose that the light degrees of freedom carry a 
small momentum q^, we can define the heavy quark momentum as: 

p^' = P^' - q^' = mqv^" + kf" , (4) 

where the "residual" momentum k^^ is defined as follows: 

fc^ = {rriH - mq)v^' - qf". (5) 
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Let us now consider the hadron scattering from an external potential. After the interaction 
we have a new hadron state, containing the same heavy quark, and carrying momentum: 



p/' = mqv''' + s*", (6) 

where is the finite momentum exchanged with the external potential in the mQ oo limit. 
If is a finite amount of momentum, then = v'^ (a finite momentum exchange cannot 
produce an infinite hadron momentum difference). This means that the velocity is a conserved 
constant of motion, i.e., it isn't any more a dynamical degree of freedom. 

When the heavy quark interacts with light degrees of freedom, only fluctuations of the 
residual momentum (of order Aqcd) have to be considered, while variations of velocity are 
certainly excluded. QCD interactions do not vary v; only weak (or electromagnetic) interactions 
can annihilate the heavy quark and create a new one that can be different in flavor, spin and 
velocity. 

In the Georgi's paper it is therefore introduced a supers election rule for the velocity of an 
heavy quark: the flelds describing the heavy quark in the eflFective theory should be hv{x) flelds, 
i.e., depending on x and v. For different v''s we have different heavy quark flelds. 

We have discussed the flavor-spin symmetry adopting the hypothesis of an heavy quark 
at rest. But now we can observe that the flavor-spin symmetry connects two heavy hadrons 
containing different heavy quark flavors only if they have the same velocity. In other words the 
SU{2Nh) flavor-spin symmetry, where Nh is the number of heavy flavors, transforms meson 
states Mq. in meson states Mq., having different i ^ j flavors, provided that Qi and Qj have 
the same velocity (not the same momentum: therefore flavor-spin symmetry is a symmetry of 
certain matrix elements, not an S-Matrix symmetry). 

Importantly, flavor-spin symmetry is not an exact symmetry since the heavy masses aren't 
infinite: the technology allowing to compute the l/mg corrections is the HQET (Heavy-Quark- 
Effective-Theory) . 

In this paper we will make frequent use of the heavy quark effective propagator. This is 
derived by the QCD fermion propagator adopting the momentum formula introduced above: 

= rnQV^" + k^". (7) 

In the niQ oo limit, the propagator: 

.7 -PQ + ruQ 



1 1 ' 



(8) 



becomes: 



2 vk' 

where we have used the relation k ~ Aqcd- The vertex describing the heavy quark-gluon 
interaction in QCD is: 

-i91>.T\ (10) 



where T"" is a generator of SU{?>)c and g is the coupling constant of strong interactions. Due to 

2 



the structure of the propagator, the vertex is always intermediate between the ^"^J'^ projectors 



and therefore the vertex in the effective theory is: 



-i.i^7.iip;T« = -i,.„T«l±l:^. (11) 
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The iij^ projectors, that appear in propagators and vertices, can be brought on the external 
Unes of Feynman graphs, where they are annihilated by the on shell heavy quark spinors 



having the property that 7 • vhv = hy, see |2.2.1| . We have therefore obtained the following two 
Feynman rules: 

i 

propagator = (12) 

V ■ k 

vertex = -igv^'T". (13) 

Since here the heavy quark mass is absent, the flavor symmetry is manifested. Moreover there 
are no Dirac matrices, therefore also the spin symmetry is manifested. 

2.2.1 l/"iQ expansion 

The Feynman rules given above can be considered as the basic definitions of the heavy quark 
effective theory. The same results can also be obtained using a field theoretical approach [ p^ ] 
avoiding the use of QCD Feynman rules like the propagator expression (^. Among effective 
theories, HQET has a particular role: one of the main goals of HQET is that of describing 
the properties of heavy hadron decays, therefore, even if there is the large separation of scales 
above mentioned, we cannot remove completely the heavy quark state from the low energy 



effective theory, see section 2.1. What can instead be eliminated in the effective theory, are the 
components of the heavy quark spinor describing its fluctuation around the mass shell since, 
in the mq 00 limit, the heavy quark is almost on shell and carries almost the entire hadron 
momentum. It is therefore useful to decompose the heavy quark QCD spinor, Q{x), in its small 
and large components: 

Q(x) = e-'^^Q^'^iH.ix) + K{x)), (14) 

where: 

H,{x) = e^-Q--l±^Q(x) (15) 
h,{x) = e^-Q--iZ-LZQ(2;). (16) 

Two properties are evident: 7 • vhy = and 7 • vHy = —Hy. Moreover, reminding the 70 
structure, one can see that in the rest frame, v = (1,0,0,0), corresponds to the upper 
components, the so called large components of the quadrispinor Q, while correspond to its 
inferior components, the so called small ones, hy annihilates an heavy quark having velocity v, 
Hy creates an heavy antiquark having velocity v. Let's consider the following virtual process 
discussed by Neubert |l^: an heavy quark propagating forward in time, at the event a inverts 
his way in the opposite temporal direction and from the h event on, it propagates again forward 
in time. In a we have the annihilation of a virtual heavy quark-antiquark pair and in b the 
creation. The energy in the intermediate virtual state, the one propagating between a and 6, 
is certainly larger, with respect to the initial one, of about 2mQ. 

Therefore, this intermediate state can only propagate over distances of order l/2mQ, which 
are very short if compared with the typical hadron physics distances, of order 1/Aqcd- The 
intermediate virtual state, that is evidently connected to the action of Hy in 6, can be simply 
substituted by the propagator l/2mQ. We can therefore state that there is no sufficient energy 
to create a virtual heavy quark-antiquark pair in HQET or, in other words, this process is 
suppressed at order l/mg. We must therefore proceed to the systematic elimination of the Hy 
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field from the effective Lagrangian therefore obtaining a non-local effective action for h^. This 
can be expanded in a l/mq series of local operators. 

In terms of h and the QCD Lagrangian for heavy quarks takes the form [16|: 

>Cqcd = Q(«7 • D - mQ)Q 

= hy{iv ■ D)hv — H^(iv ■ D + 2mQ)H^ 

+ K{ij-D)H^ + H^{ij-D)K, (17) 

where: 

jji' = D^" - yf'v ■ D. (18) 

We conclude that h describes massless light degrees of freedom, while H describes fluctuations 
having a mass of 2mQ. The latter must be eliminated from the effective theory. From the last 
two terms in >CqcDi describing the creation (annihilation) of quark-antiquark pairs, we can see 
that H and h fields are mixed together. If we compute the functional derivative of £qcd with 
respect to H, we obtain the following equation of motion for H: 

{ivD + 2mQ)H^ = i-fDhy, (19) 

that can be formally solved for and the resulting expression can be inserted in £qcD; giving: 

-Ccff = K{iv ■ D)hv + Kiil ■ D- -17 • Z))/i„, (20) 

2mQ + IV ■ D 

where the second term describes the virtual process discussed above. In the momentum space, 
derivatives acting on h fields correspond to powers of the residual momentum A;, therefore we 
can perform the following power expansion: 

1 °° / iv 
CeS = K{iv ■ D)hy + "^Kij-Dl-- j ■ Dh^. (21) 

It is not difficult to prove the following identity: 

1 + 7 • f 



(22) 



where G^u is the gluon strenght tensor field and the well known property [iDfj_,iD^] = igGfj_iy 
holds 1 26]. Considering the n = term in the expansion (pl|), one finds the interesting result: 

£eff = hviiv ■ D)K + -^K{iDfK + -^Ka^^G^^K + 0{l/ml). (23) 
2mQ AniQ ^ 

The mq oo limit selects only the first term in the preceding Lagrangian. Since the heavy 
quark mass is large, but not infinite, all other terms are to be considered as corrections, that, 
as we can see, are included in the effective theory in a systematic way. The first term in Ces 
allows to write down the Feynman rules for propagator and gluon vertex already discussed 
before. Let us rewrite it including a sum over heavy fiavors and a sum over velocities: 

4ff = E/ d'^r—^\iv ■ D)h^:\ (24) 

i=l 
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Since in this Lagrangian there aren't terms containing mg, we can deduce that is invariant 
under flavor space rotations. Moreover, since no Dirac's 7 are present, the interactions among 
heavy quarks and gluons conserve heavy quark spins. This is the SU {2Nh) symmetry aheady 
discussed. To be rigorous, since Lorentz transformation can boost the heavy quark velocity, 
the symmetry group should be Lorentz x SU{2Nh)°° ■ 



Let's now consider the two operators of order l/mq in (24). To understand their role it is 



convenient to write them in the rest frame of the heavy quark v = (1, 0, 0, 0): 

^ -K{ibfK ^—K{iiyfh^ (25) 



2mQ 2mQ 

^ -Ka^.G^^'^K ^ -^K{S-Yl,)K, (26) 



Amq mg 

where H* = —^e^^^G^^ and S is a spin operator defined as a 4 x 4 matrix with Puli matrices 
^ on the diagonal. Therefore the first operator is a kinetic energy operator connected to the 
residual motion of the off-mass-shell heavy quark. The second operator is the non Abelian 
extension of the Pauli interaction describing the chromomagnetic heavy quark spin coupling 
with the gluon field. We find confirmation that the heavy quark spin is decoupled by a factor 
of l/mq. 

2.2.2 Relations with QCD 

To match HQET with QCD at high energy, one must include some corrective effects in HQET 
due to high energy virtual processes. For example, the weak current transforming the flavor 
from 6 to c must be corrected at the Og order both in QCD and in HQET. We can expect that 
there are differences between these corrections. These differences instruct on how one should 
modify the coefficient of the weak current in HQET and on what terms should be added to the 
HQET current to guarantee the correct matching between the low energy and the fundamental 
theory. Let us consider for example the case of the current b-y^c. The 'j^ of QCD has to be 



substituted by the of HQET where pi: 



F^ = (1 + Co^) r + ^Yl ^^^r , (27) 

i 

and F^ are new structures containing v and v' , i.e., the velocities of the heavy quark before and 
after the weak interaction vertex. In practice, this type of calculation is made by comparing 
the vertex diagrams where the fermionic heavy quark current is coupled to the weak current, 
up to order q^. We have therefore to compare four QCD diagrams with four HQET diagrams. 
The difference between the two set of diagrams lays in the Feynman rules describing the strong 
vertices and the heavy quark propagators. The first of the mentioned four diagrams is the three 
level diagram (the simple tree weak vertex). In the remaining three diagrams, one should close 
the gluon line on the heavy quark line according to the three diagrammatically possible ways. 

2.3 Chiral lagrangians 

The QCD Lagrangian with three massless quark flavors incorporates the U{3)l x ^^(3)_r global 
symmetry [^. The left- and right- handed components ql = {ul, (11, sl) and qr = {ur, dji, sr) 
transform respectively as the fundamental representations of U{3)l and U{3)r respectively. 
Anyway, the symmetry group of the quantum theory is a subgroup of U{3)l x U{3)r, this is 
what is usually called anomalous breaking at the quantum level of a symmetry of the classical 
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Lagrangian. The reason for this phenomenon is that U{1)a is not a good symmetry of the theory 
since its generator, Q5, is not a time independent quantity due to the presence of instantonic 
configurations of gauge fields in Yang- Mills theories. The quantum theory has therefore the 
following symmetry group: SU{3)l x SU{3)r x U{1)v, the chiral symmetry. Anyway the 
physical states are invariant only under SU{3)v x U for example the baryon spectrum is not 
doubled in two spectra having opposite parity, but it is well described as an octet representation 
of SU{3)v having baryon number 1: this is the phenomenon of spontaneous symmetry breaking. 
The mechanism underlying the spontaneous symmetry breaking of chiral symmetry is most 
likely of non perturbative nature. The energy scale associated to this phenomenon is ~ 1 
GeV (this point will be discussed with greater detail later). 

Each broken global symmetry implies the exitence of a Nambu-Goldstone (NG) boson 
emerging as a scalar massless particle induced by the non symmetric structure of theory's 
vacua. The chiral symmetry is not an exact symmetry because light quark masses are not 
exactly zero: they are simply small if compared with Aqcd- Therfore we should expect to have 
pseudo-NG-bosons having small masses. 

Light quark masses are slightly dissimilar to each other causing that NG bosons will also 
have slightly dissimilar masses. On the other hand, since light quark masses are much smaller 
than AqcD) the light baryons are almost all degenerate in mass (because, differently from NG 
bosons, sending to zero the light quark masses, baryon masses should tend to a value different 
from zero). 

Due to the lightness of pseudo NG bosons, it emerges a hierarchy of energy scales allowing 
to decide that NG bosons interactions at energies much lower than A^ can be described within 
a chiral effective theory. Even if in this case the separation among energy scales is not as large 
as in the case of heavy quark effective theory, the chiral effective theory, developed in seminal 



papers by Weinberg ||lj], Manohar and Georgi |2^] and by Gasser and Leutwyler |3^, is a great 
success. 

Since chiral symmetry is spontaneously broken, we have a chiral condensate different from 
zero that we can write as pl| : 

{i^jRi^iL) = cT^ij, (28) 

where c is the value of the condensate, while S defines a direction of the condensate in the flavor 
space. All S's orientations are degenerate vacua that are mapped one in the other through the 
SU{3)l X SU{3)r transformations: 

E ^ LSi?+. (29) 

S is normalized in such a way that S^S = 1. 

When S has a spatial dependency, then we are dealing with a Goldstone excitation: low 
energy excitation are in fact characterized by a vacuum configuration that varies from point 
to point in the space, being the orientation of the vacuum state a function having a smooth 
dependence on the position (think of spin waves in a ferromagnet). The excitation energy is as 
small as one likes when one considers very small S(x) variations over large lenght scales: this 
is the case of NG bosons. 

In order to construct a chiral effective theory, one needs to follow the instructions that 
we have already described for a general effective theory. One must write the most general 
Lagrangian, containing the T, field, consistent with relativistic invariance, VCT, QCD chiral 
symmetry. 



In the chiral effective Lagrangian there aren't terms not containing derivatives |13]. If there 



was such a term, it could only be constant (think of Tr[S+S] where S+S = !)• Every invariant 
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function of S without derivatives, is just a constant. If we don't exclude terms not containing 
derivatives, we could have a Lagrangian containing only zero momentum S fields. But a NG 
boson with zero momentum simply does not exist, it is equivalent to the vacuum. 

The first non trivial term that we are going to consider in the chiral Lagrangian construction 
is the one having two derivatives: 

£2 = ^Tr[a/,E9'^E+], (30) 

this is the only allowed one with two derivatives, /tt is the pion decay constant, i.e., one of the 
energy scales characterizing the pion (the other one is its mass). The E field is described as 
the exponential of NG boson fields: 

S(x) = e~S^, (31) 

where: 

Ti{x) = 7r'^(x)r^ (32) 

r° being the 8 generators of S'C/(3)-colour. Such an expression for the S field amounts to 
consider the NG fields as the angular variables describing the vacuum rotations. Moreover let's 
observe that under the transformation (^), the vr fields transform non linearly as complicated 
functions of L and R. The exponential representation (^), is only a particular one, the 
simplest, among all the possible non linear representations of tt fields. All of those give the 



same S-matrix |32] 



Let us now take the first term introduced in (|30|). This can be expanded in powers of the 
pion momentum in the following way: 

C2 = Tr[9^7ra'^7r] + ^Tr[[^, d^7:f] + , (33) 

where higher order terms give non-linear interactions with NG bosons. Besides C2 we can 
introduce higher dimensional operators Cn, i-e., containing a larger number of derivatives. 
Manohar and Georgi have shown that the energy scale A controlling this expansion in an 
increasing number of derivatives ^ is A = A-^ and they estimate A^ ~ 'i'T^fn — 1 GeV. This 
limits the range of validity of the Lagrangian to the low energy domain since the pion momenta 
have to be small compared to A^ (otherwise the derivative expansion gets divergent very soon) . 
We will make these points clear in the next sub-session. 

2.3.1 A^ 

As stated above, besides (|30|), the chiral effective Lagrangian may also contain terms with an 
higher number of derivatives. Let's consider £4 containing four derivatives, i.e., proportional 
to: 

£4=^Tr[9^Sa^S+a''Sa^S+]. (34) 

An higher number of derivatives means an higher number of pion momenta. As already stressed, 
in a low energy effective theory a term such as £4 needs to be multiplied by some inverse 
power of A, an energy scale characterizing the upper energy bound of the effective theory and 
representing the parameter that controls the convergence of the pion momentum expansion. 
This A can be associated to the energy scale of spontaneous chiral symmetry breaking. Ay = A, 
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and this association will allow for an estimation of that is going to be used throughout this 
work. We will follow an argument due to Manohar and Georgi [^9|. We can normalize £4 
relatively to the lowest order term (|30|), adding two powers for each additional derivative. 
The term having the correct mass dimension is: 

£4 = ^TV[9;,S9^S+9''Ers+]. (35) 

X 

If was much higher than pseudoscalar masses, then all orders having an higher number of 
derivatives with respect to (|30|) would be negligible. But we cannot formulate this hypothesis 



since radiative corrections to (^OD produce a term as (|35D and also higher order terms having 
an infinite coefficient. Thus, even if these terms are absent or negligible for a certain choice of 
the renormalization scale, they can be present and important for a different one. On the other 
hand, one could reasonably think that the /^/A^ coefficient should be numerically larger or 
equal to the variation induced in it by an 0(1) shift in the renormalization scale of radiative 
corrections to (pO|). 

These considerations get particularly clear if one examines the vr — vr scattering process. If 
we refer back to (|30|), we see that each four pion vertex has the following structure: 

J IT 

Using two of these four pion vertices to generate the one-loop diagram shown in fig. ^, we 
can obtain the particular case of pion-pion scattering diagram with all derivatives acting on 
the external legs (the same diagram in fig. ^ may have only two derivatives acting on external 
legs, while the two remaining could act on the internal ones. In such a situation we are facing 
a quadratically divergent diagram). 

If instead one refers directly to (|35[), the fundamental four pion vertex has the form: 

4 4 

f2A2' y^'l 

therefore one can also write the pion-pion scattering diagram with all derivatives acting on 
the external legs simply using the tree level diagram extracted from £4. The diagram in fig. 
^, generated by two insertions of C2, is just a one-loop correction to the tree level process 
described by a single insertion of £4. The diagram in fig. ^clearly gives: 



dividing by the correct symmetry factor and introducing a cutoff Aco and the renormalization 
scale fi. To be consistent with a chiral effective theory scheme, we should limit the momenta 
circulating in the loop to Aco — A^^. Let's now rescale /i in (^8|), of an 0(1) quantity: let's take 
the Neper number. We obtain a variation of (p8t) amounting to: 



1 p^TT^ , , 

(39) 



(4vr)2 /4 



The shift in the renormalization point ^ can be absorbed into a redefinition of the coefficient 
/^/A^ in (35) and we have that ( |39| ) corresponds to a change in that coefficient amounting to: 
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Figure 2: One loop diagram for vr — tt scattering. 



But, as above observed, we can expect that: 

f2 1 

— > (41) 

therefore: 

Ax < ^ttU, (42) 

which suggest to use, as an estimate of the energy scale associated to spontaneous chiral 
symmetry breaking, the following one: 

= 4ttU ^ IGeV. (43) 

In this work we will make use of a cutoff value close to A^ for the same physical motivations 
here described. The fact that /tt is a measure of the strength of the symmetry breaking is also 



discussed in [33|. 



2.3.2 The Manohar-Georgi Lagrangian 

The effective Lagrangian defined below the chiral symmetry breaking scale contains, besides 
pion fields, also light quarks and gluons. Let's define the ^ field in the following way: 

C = e'^, (44) 
15 



i.e., ^ = vl]- We know that under SU{3) x SU{3) transformations of the S field, the NG 
bosons, represented by the vr matrix, undergo a non-hnear transformation where tTq — > vr^, vr^ 
being non-linear functions of tt, L, and R. The transformation properties of vr fields determine 
also the transformation properties under SU{3) x SU{3) of the ^-fields. Since S = 

C^^' = L^U+ = U^R+, (45) 

where [/ is a non linear function of L, R and vr that can be written as an ordinary SU (3) 
matrix in the following way: 

U = e'\ (46) 

and the Hermitian matrix v contains "non-linearly" the SU{3) x SU{3) symmetry. The U 
matrix is invariant under parity transformations exchanging L with R and vr with — vr. If 
L = R, the chiral transformation reduces to a simple SU{3) transformation and U = L = R. 
Through ^ matrices we can construct two auxiliary fields: 

= h^+d'^^ + ^d>^e) = ^[^,d^7r] + (47) 



A'' = -ktd^C-^d'^t) = ^d,7r + , (48) 



which, under chiral transformation, behave like this: 

V^* ^ UV^U^ + Ud^U+ (49) 
^ (50) 

In particular, due to the transformation property of V^, we can treat it as a gauge field in a 
covariant derivative: 

D^' = ^^' + V^'. (51) 

Let us now form the Lagrangian terms related to n, d and s quarks considered together in 
a unique triplet V of flavor-S'C/(3) . V is supposed to transform under SU{3) x SU{3) in the 
following way: 



The only term without derivatives is: 



ip Uip. (52) 



-rm/jip, (53) 



where m is not the current mass of the QCD Lagrangian. m is a constituent mass whose origin 
can be related to the spontaneous breaking of chiral symmetry (we will come back on this point 
when we will introduce the constituent light quark mass in the CQM model). 
We also have two derivative terms. The first one is the kinetic term: 

iip'j ■ Di/j, (54) 

the other one is the interaction term between the light quarks and the pion. We will use the 



PCAC language |34] to introduce this term. The derivative operator of the axial current, dA, is 
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(a) 





7Z 








(b) 



Figure 3: (b) represents the pion pole contribution to the diagram in (a) 



a pseudoscalar operator with odd G-parity, isospin one and hypercharge zero. It has therefore 
all quantum numbers of the i-component of pion triplet. Let us consider the matrix element: 



(55) 



where ■0 is a light quark field, see fig. ^a). Many diagrams contribute to diagram in fig. |3|(a); 
we choose that in fig. ^b). This is equivalent to: 



q2 _ ^2 



where Amp{^ ip'-rr) denotes the amplitude for the process ip — > ip'n (the i comes from the 
propagator). Clearly: 



Amp(V' ^ V'vr) = ^-^{i^'\dA\il;), 



(57) 



in the — > limit. PCAC hypothesis consists essentially in defining the following off-mass- 
shell amplitude: 



(58) 
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where q is the pion momentum. At this point one must suppose that this off-mass shell 
amplitude varies smoothly within the (0,m^) g^-range. The chiral limit of Amp is defined 
taking the m-,^ — > limit and then the mass-shell-limit — > (the two procedures do not 
commute) hoping that the hypothesis on the smooth variability with respect to holds well 
(as is confirmed by the Goldberger-Treiman relation and by its physical consequences). The 
chiral off-mass shell amplitude is therefore: 

A^(V' ^ V^'tt) ^ -^{ip'\dA\i:) = ^(V^'I^'^IV'). (59) 

Jit Jtt 

We can therefore have a Yukawa type coupling in the Lagrangian: 

^ _ 

-^V'T^TsV^^'^TT, (60) 

Jtt 



and generalize it, see (48), in the following way: 

tpj ■ A'y^ip. (61) 

The Manhoar-Georgi Lagrangian to lowest oder, including also colour in the covariant deriva- 
tive, is: 

^eff = i'i'j ■ (iD + V)]-ip + gAi'l ■ Aj^'ip - rri'ilt'il; 

+ :|TV[a''Ea^S+] - iTV[G^,Gn, (62) 

where G^^ is the well known non-abelian strenght tensor of the gluon field = G^T"". qa must 
be computed through a matching with QCD or extracted by a comparison with experimental 
data, or, as will be our case, with some more fundamental effective model. 

2.3.3 Heavy mesons and chiral Lagrangians 

We are going to discuss of an effective chiral Lagrangian describing the interaction of soft pions 
(and kaons) with mesons containing an heavy quark. As we saw in section |2.2| , heavy meson 
fields should be described through the HQET formalism. In order to implement the flavor-spin 
symmetry, the field describing an heavy meson has to be independent on heavy quark mass 
and spin. On the contrary, it can be characterized by the total angular momentum of light 
degrees of freedom. To each si value, there corresponds a doublet of states degenerate in mass 
with total angular momentum J = it ^. In correspondence of = ^, for example, we have 
the P and P* mesons being respectively the pseudoscalar and vector components of the spin 
symmetry doublet. If the heavy quark is c, P and P* correspond to D and D* , while if the 
heavy quark is 6, they are the states B and B* respectively. 

Let us consider the negative parity doublet (P, P* ) . We can associate a unique super-field 



H [35| describing both states. This super-field must have two spinor indices: one connected to 
the heavy quark and the other to the light quark. The structure of H is that of a 4 x 4 Dirac 
matrix. If one performs a Lorentz transformation, H behaves like: 

H ^ D{A)HD{A)-^, (63) 

where A is the usual 4 x 4 representation of the Lorentz group. An explicit matrix representation 
is the following: 

H = ^^:^[p:r-P75], (64) 
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and we define: 



H = 7o^^o. (65) 

V is tlie velocity of the heavy meson and the following transversality condition holds: v^P* = 
and Mh = Mp = Mp*. We mention also the following useful relations: 7 • vH = —H'j ■ v, 
H'f ■ V = —7 • vH = H. P and P*^ are the annihilation operators normalized in the following 
way: 

(VAC|P|Qg(0-)) = ^fmi (66) 
(VAC|P*'^|Qg(l-)) = e^v^- (67) 

The formalism apt to describe higher spin meson states and the effective Lagrangian terms 



associated to them, is extensively developed in |3G]. We are interested in considering the p- 
wave {i = 1) of the Qq system. HQET predicts the existence of two degenerate doublets: 
(O"*", l"*") and (l"*", 2"*") for each heavy quark c or b. The related superfields are: 



1 + 7 • f 

T 
1 + 7-7; 



S = —i—[Pi',ri6-Po] (68) 



(69) 



These two doublets have respectively se = ^ and = |. This classification with respect to Si 
is the more reasonable one since we know that the dynamics is completely independent on the 
spin and on the mass of the heavy quark Q. Observe that in the limit of infinite heavy quark 
mass, S£ and sq are separately conserved. We can introduce the total spin J, defined as the 
total angular momentum of the heavy and light quark in the rest frame of the heavy quark: 

J = + SQ. (70) 

Heavy mesons can interact with vr fields through their light degrees of freedom. The inter- 
action terms of the NG boson octet with heavy meson fields must be written in an effective 
Lagrangian including the essential symmetry properties: chiral symmetry and heavy flavor-spin 
symmetry. This heavy-light chiral effective Lagrangian can be expanded with respect to: 

• NG bosons momenta 




powers 



An heavy-light Lagrangian has been introduced almost simultaneously by different groups 

C = ^TV[9^S9^S+] - Tr[Hiv ■ DH] + gTj:[HH-f ■ ^75] + , (71) 

where ellipses indicate the presence of an infinite number of operators having higher dimension- 
ality, including those responsible for explicit chiral symmetry breaking, i.e., terms containing 
light hadron masses, and those of order Aqcd/w-Qi violating the flavor-spin symmetry (such 
as the color magnetic moment operator). The covariant derivative has been defined in dST]). 
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Let us consider for example the term having the factor of g in eq. (|7l|): it describes the 
coupling of i^-type mesons with NG bosons, see eq. (^8|). We will come back on terms of this 
kind many times. Including also S and T mesons we have: 

C = igTT[HH^ ■ ^75] + ig'TT[SSj ■ ^75] + ig"Tr[f^'T^-f • ^75]. (72) 

At the lower order of the derivative expansion we can also write the interaction terms describing 
transitions between different doublets; for example: 

C = ifTT[ST^'A^-f5] + if'TT[HS-f • ^75] + h.c. (73) 

A particular case is that of transitions between T and H states. Let's suppose to consider the 
following s-wave interaction Lagrangian: 

£ = irTJ:[HT^'A^-f5]■ (74) 

We can therefore write the S-matrix element: 

(outL>7r|L>2m), (75) 



as: 

{D-K\iC\D2) = -r(L>|P|VAC)(7r|^|VAC)(VAC|P2^n^2)Tr 

= iry/mDmD^^2e^''vu, (77) 

where the first order in the expansion ( ^8|) has been used. Due to transversality, this term is 
certainly zero, since e^^'^Vy = 0. Since the process we are discussing is entirely due to the strong 
interaction, the velocity conservation rule, introduced in |2.2| , holds, i.e., the velocity of D and 
D2 are the same. 

We conclude that the Lagrangian ( [7^ cannot be the right one to describe the T Htt 
transition. We need one more Lorentz index coming from the insertion of another derivative 
under the trace sign. This derivative should be accompained with a negative power of A^, 
giving the right mass dimension to the interaction term and controlling the expansion in NG 
bosons momenta. 

The d-wave Lagrangian is: 

C = ^TTlHTf^iiD^j ■ A)j5] + ^Tr[^r^(i7 • DA^H] + h.c. (78) 
With an analogous approach super-fields having higher spins may be constructed. 



1 + 7 • 7; 



-7^ 



(76) 



3 CQM 

3.1 The CQM model 

CQM is a Constituent-Quark-Meson- Model that has been introduced and discussed in a number 
of recent research papers [H, H, H, ||. The model, based on an effective Lagrangian 
describing quark-meson interactions, is relativistic, incorporates the flavor-spin symmetry in 
the heavy sector and the chiral symmetry in the light quark sector. In the following sections 
CQM will be reviewed in detail. Section 4A. and subsequent sections are instead devoted to 
the CQM phenomenological applications to heavy meson physics. 
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In the weU known old paper [43|, Y. Nambu and G. Jona-Lasinio discuss the possibiUty 
that the nucleon mass can be due to an unknown primary interaction bounding hypothetical 
massless primary fermions. In their model, the same interaction bounds nucleon pairs giving 
rise to pions and the mass of the Dirac particle emerges as a result of the primary interaction 



in the same way as the energy gap in a BCS superconductor |44] is connected to the formation 
of correlated Cooper pairs of electrons, as a result of a phonon-mediated "primary" interaction 
(the finite energy that is needed to break a Cooper pair is proportional to the BCS gap). 

The primary interaction in the Nambu-Jona-Lasinio model is a non linear four fermion 
interaction, already discussed in an older paper by Heisenberg et al.: 



(79) 



In some recent papers |45, the problem of the bosonization of a Nambu-Jona-Lasinio 
(NJL) like Lagrangian, where the fundamental fields are quarks, has been extensively studied. 
The aim is to derive an effective theory of mesons starting from a model (the NJL model) 
incorporating global chiral symmetry and its spontaneous breakdown: order parameter of the 
chiral symmetry breaking is the light constituent quark mass which is proportional to the chiral 
condensate and can be calculated through a gap equation [^], as we will briefly see. 

Let us consider a U{Nf) multiplet '0 = V'/,c) the A matrices A G U{Nf), normalized accord- 
ing to Tr/A°A^ = 26°^^, and let's write: 



N--1 

E 

a=0 



- A" 



+ 



(80) 



(Interactions (|79|,80) are interconnected by Fierz theorem). Taking the derivative of (| 
respect to -0, we can write the following equation of motion: 



with 



G 



(81) 



where we require that (tpX^'i^^ip) = [47|. If we also require that [iJjX'^iJj) = for each a ^ 0, 
following an analogy with the approximations made to solve the BCS equation of motion 



44], and remind that A" = ^/-^l, defining: 



^(VACIV^VIVAC) = -mdyn, 



(82) 



the equation of motion becomes: 



{i-f -d- mdyn)V' = 0. 



(83) 



Let us now consider that: 



G 



N 



f 



(VACIV^VIVAC) 



G 



lim TV/Trc(VAC|V'/,c(a;)?/'/,c(0)|VAC) = iGNcApia 



(84) 



where Af{x) is the Feynman propagator: 



7-p + mdyn 



(85) 
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Using the Dirac equation, we obtain the fohowing expression for the dynamical mass: 



-dyn = r d^p ^ "^"r ^ . . (86) 

(27r)4 J pi - m^y^ + le 

The non-trivial solution, mdyn 7^ 0, is connected to the spontaneous breaking of chiral symme- 
try. 



In the papers by D. Ebert et al. [45], the NJL model is extended to include two light quark 
fields and an heavy one. In such a way the bosonization produces collective meson fields having 
a light-light or heavy-light constituent quark content. 

3.1.1 Bosonization 

The basic idea of the bosonization technique is that of re-formulating a field theory written in 
terms of microscopic degrees of freedom, such as quarks and gluons, as a field theory in which 
meson fields are on the same footing of elementary fields. Many attempts to bosonize the QCD 
Lagrangian, i.e., to yield a meson theory mathematically derived from first principles, have been 
unsuccessfully performed. Some progress in this direction has been made in two-dimensional 
QCD |§. 

In the path integral language, this is how bosonization works: 

j DqDqe' f-^'^'^ j DaDTTDp...e'^-^^°% (87) 

where Da, Dtt, Dp... are the integration measures associated to the meson fields. The effective 
Lagrangian £bos is written as a function of these fields. 

The Hubbard-Stratonovich transform is the first step in (|87|): the four quark NJL interaction 
is substituted by Yukawa couplings of the quark fields with meson fields: 

e^L-CNjL^ I DaDi:...e'-^^^'^3-L. 



£'p^jL is the semi-bosonized £njl- ^njl Gaussian with respect to functional integration 
over microscopic fields. Therefore, integrating over Dq and Dq one obtains a determinant 
containing the meson fields. This can be loop expanded and the Feynman diagrams coming in 
this expansion can be evaluated in the region of small meson momenta, the most interesting 
for our purposes. 

The NJL Lagrangian studied in |45] is: 



/:njl = (v^Tmy^) (v^^'^y^ ) ' (89) 



where q = {u,d,s)'^, Qy = b or = c, tp = {q,Q)'^ and G is a coupling having dimension of 
(mass)"'^, while A are matrices of SU{Nc). The free Lagrangian is therefore the sum of the 
familiar Dirac Lagrangian for light quarks and the free effective Lagrangian for heavy quarks: 

Co = q{i^ ■ d - 'ih)q + Qy{iv ■ d)Qy. (90) 

The bosonization is then performed on: 

£ = £o + >Cnjl. (91) 
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Fierz theorem allows to rearrange the bosonized Lag rangian iZ^os ^ sum of three pieces. We 
are interested only in two of them: £bos = + -C'*'- The former is related only to the light 
degrees of freedom, the latter includes also the heavy quark and heavy meson fields. The third 
term, C^^, is not relevant when one is interested in studying the physics of mesons containing 
only one heavy quark, as is the case here. 

£ in eq. (plj), has the global colour symmetry, the chiral SU{3) x SU{3) symmetry (as the 
current light quark masses go to zero) and the flavor-spin symmetry of HQET (observe that 
also the interaction term is independent on the heavy quark mass and spin). These properties 
are preserved in £bos- 

Technical details about the bosonization method are far beyond the scopes of this work. 
The interested reader is referred to ||4^ and references therein. 
The CQM La grangian is a phenomenological extension of jC^qs- 



3.1.2 The CQM efTective Lagrangian 



As discussed in the last section, the CQM Lagrangian is made up of two terms, like £bos) but 
does not exactly coincide with it (£cQM 7^ ^bos) for reasons that will be clear soon: 



CQM 



hi 



(92) 



The first term describes the light degrees of freedom and it is very similar to the Georgi- 
Manohar Lagrangian given in (|6^). The differences are that, in the CQM Lagrangian, there 
are no gluons and the light fields are defined differently: 



= xh ■ (id + V)]x + XT • ^75X " mxx + ^TV[a'^Ea^E+] 



(93) 



where now we define /tt = 130 MeV. The absence of gluons is rather plausible since this 
Lagrangian originates from the bosonization of an underlying NJL interaction Lagrangian, 
where gluons are absent from the start. 

The light x fields are also a consequence of the bosonization of an underlying NJL. What 
emerges is that x = where q are the familiar light quark fields and ^ = e f . We will always 
consider the expansion of ^ to be truncated at the zero-th order in the pion field. We will need 



the first order of this expansion only in section 5.4 



From a detailed comparison with (|62[), it is also evident that in (93) = 1, again as a result 
of bosonization. The mass m in (p^) is dynamically generated according to the mechanism 



explained in section 3.1 



Let us now focus on C^K Here we have the Yukawa type interactions between quark and 
meson fields emerging from bosonization, plus two phenomenological terms put by hand: 



+ 



x[H + S + if^'^ ] Q^ + h.c. 



^ TT[HH] + ^Tr[SS] + ^TT[f^^T,]. 



2G3 



2G4 



(94) 



The first term is the well known heavy quark kinetic term of HQET, see 2.2 

The second term is responsible for the Q — Meson — q vertices shown in fig. ^j: these are the 
most relevant aspect of the CQM model. The meson fields H, S and T have been introduced 
in 



2.3.3. 



The vertices for H and S mesons have been derived from bosonization. The vertex involving 
the T field is instead a phenomenological term, introduced according to the philosophy of 
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Figure 4: The CQM Q - Meson - q vertex. 



effective theories (A = 1 GeV in eq. (p4D). This value of A should also be assumed as the 
ultaviolet cutoff in the regularization of the model A 



A;^-, see the discussion in |2.3 .1 



Following [45| we will adopt as UV cutoff A = 1.25 GeV. This is a value quite close to the A^ 



discussed in 2.3.1. On the other hand one finds that the CQM phenomenological predictions 
are not very sensitive to the variation of the UV cutoff, at least if one varies it within 10— 15%. 
Last three terms in are the kinetic terms for H, S and T fields. The first two terms come 



from bosonization, while the third is inserted by hand as a phenomenological term resembling 



the first two. Bosonization also predicts that 2^ 



'2GI, 



and we can therefore write: 



C 



hi 



+ 



Qv{iv ■ d)Qy 
1 



X\H + S 



- 9^ 

1 



h.c. 



2G. 



Tr[(i/ + 5)(/7-5)] + — TV[r^r, 



(95) 



is crucial for the CQM calculation of the coupling 



The dynamical information 2^ = 
constants: there are not sufficient experimental data to determine two different constants G3 
and Go. 



The kinetic terms will be rewritten, in the form discussed in [35|, in the next section, where 
we will also discuss the problem of determining the mass difference between the H and S 
multiplets. Again the dynamical information 2^ = 



^ — helps this evaluation. We will call 
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Ah, As and the mass differences between tfie masses of H, S and T multiplets and the 
heavy quark there contained. The mass difference between S and H will be simply A5 — Ah 
and it will be zero as soon as the light constituent quark mass m — > (i.e., in the chiral 
unbroken phase). 



3.1.3 Regular izat ion 



As we pointed out in |3.1.2| , CQM is the fusion of a Manohar-Georgi like Lagrangian for the 
light quark sector with a quark-meson Lagrangian for the heavy quark sector. We therefore 
know that the upper energy scale, i.e., the energy scale over which the effective theory should 
be substituted by a more fundamental theory, is A^. It could seem strange that the heavy quark 
mass is itself higher than the UV cutoff, but we have to remind that in HQET the on shell 
momentum of the heavy quark, mQV, is not a dynamical quantity since, due to the velocity 
superselection rule, v is not dynamical, see |2.2| . The dynamical quantity due to the interaction 
between the heavy quark and the light degrees of freedom is the residual momentum k^, which 
is necessarily k ~ Aqcd < A^^. 

CQM does not include the gluon fields, as it is obvious considering that it results from a 
path integral bosonization of a NJL model, and does not incorporate confinement of quarks. 
This may appear as a strong limitation of the model but, according to a common opinion, 
it is physically much more important to work with non confining models possessing chiral 
symmetry and its spontaneous breakdown than with confining models where chiral symmetry 
and its breakdown are not properly incorporated. In the former case one is describing a world 
which is essentially the same as the real one for what concernes the hadronic spectra; the only 
difference would be that of the theoretical admissibility of asymptotic quark states. The latter 
case presents an hadronic spectrum completely messed up with respect to the observed one. 
We show here how one can face the problem that CQM is not a confining model: introducing 
an infrared cutoff fi. 

The kinematical condition for an heavy meson having mass M to decay into its free con- 
stituent quarks is: 

M >mQ + m. (96) 



Since the meson momentum is P = rngv + k, where P = Mv, eq. ( |96D is equivalent to the 
condition v ■ k > m. In the frame where the heavy meson is at rest, the latter condition means 
ko > m, i.e., inf(A;) = m. Therefore one should consider residual momenta k larger than m 
to be sure that the unphysical threshold condition ( |96|) holds, as it should in a not-confining 
model. For lower k values one is in the energy region where confinement must be necessarily 
taken into account. 

On the other hand, the value of the constituent light mass is determined by a gap equation 
(see also (||)) [||: 

(S) = m = rh + 8mli{m'^), (97) 



where the chiral ray E has been defined in (|28|), while the Ii integral is given in the Appendix 
together with other integrals met in CQM applications. Ii is calculated with an UV and an 
IR cutoff introduced according to the Shwinger's regularization method, as we will discuss in 
a while. As the infrared cutoff varies, the m value varies accordingly and, following what we 
have observed before, we can choice as an infrared cutoff /i ~ m (the running momenta in the 
CQM loops we will deal with, are of the size of the heavy quark residual momenta). In the 



second paper in reference |45|, it is shown the m vs. /i plot obtained from (|97|) for a fixed 
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A. This plot has the typical shape of a second order phase transition order parameter with a 
critical at /Xc — 550 MeV. For ^ > ^Uc, m is zero, i.e., the chiral symmetry is unbroken. For 
fi = m = 300 MeV, one is in the broken (physical) phase at the edge of a plateau. 

Therefore the boundary energy values of the effective theory are chosen to be /i = 300 
MeV, A ~ A^ ~ 1.25 GeV, and the light constituent mass is dynamically generated by a NJL 
gap equation: m = 300 MeV (which represents the degenerate u and d masses. We will not 
consider s quarks). 

The last step is the choice of the prescription to implement the cutoffs in the calculations. 
For a non renormalizable model, this step is part of the definition of the model itself. The 
proper time Shwinger regularization has shown to be the most adequate for our purposes. 

After a continuation of the light propagator in the Euclidean domain, the following pre- 
scription is used: 



d'^iE^r^ - [ d^E dse-<''E+^"\ (98) 

All the CQM calculations are performed applying this receipt. If one tries to insert the cutoffs 
as the bounds of the Euclidean integral measure, besides the problem that Euclidean translation 
invariance is then lost, one also has to face the problem that the choice of the infrared cutoff 
is not only conditioned by m (see the discussion made above), but also by A^f, that is the 
free parameter of our model. The regularization receipt (p^) acts in the sense of modifying 
the Euclidean light propagator through a factor depending on the difference of two exponential 
functions of {£\ + m^). This could affect the Ward-Takahashi relation for, e.g., the vertex of 
the axial current A with the light quarks, causing the emergence of a mass term for the pion, 
even if one considers the chiral limit dA = from the beginning. Anyway, due to the structure 
of the regularization receipt, this should be a very soft effect. 

3.2 Renormalization constants and masses 

The simplest CQM loop diagram that one can obtain contracting two vertices Q — Meson — g, 
see fig. ^, is the meson self energy diagram shown in fig. ^. 

For in and out H fields, we can write down the following loop integral: 

The rules applied are the standard ones for loop integrals. The expressions of the usual Dirac 



propagator and of the heavy quark propagator, defined in 2.2, have been inserted in the 
integral together with the vertex prescriptions derived from the heavy- light Lagrangian C^K 
The regularization procedure is that of the Shwinger's proper time, see [3.1.3 . 
First of all let us observe that we can perform the expansion: 

n(?; • /fc) ~ n(A) + n'(A)(t; • /c - A), (100) 

since we know that k smoothly fluctuates around (M — mQ)v, i.e.: 

F = Av^ - q^", (101) 

where q parameterizes this small fluctuation, see eq. (^), and A is deflned as A = M — mg 
modulo 1 /niQ corrections. This expansion of IT can now be inserted in the self energy expression 
for the H fleld (for S and T fields the procedure is exactly the same) and subtracting from 
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Figure 5: CQM meson self energy diagram. 



£ the counter-terms Tr[//n(f • — Tr[5n(u • k^S\ and — Tr[T^n(?j • A;)T^], one obtains a 

modified kinetic part of tliat can be written as follows 



£,^en = -TV[^,en(i^^-5-AH)Hren]+TV[5ren(i?^-5-A5)5rcn] 

+ TVfc(it;-9-AT)rrenM], (102) 



provided that: 



2G3 
1 

2G7 



HhCAh) = ^s{^s) (103) 
nT(Ar) (104) 



//rcn = (105) 
5ren = ^ (106) 



Tren = (107) 



where the renormalization constants Z are defined as follows 

' d_ 

dx 



z-i = f An(x) ) , (108) 



with j = H, S,T. As showed, the kinetic part of that originally was written as: 

_i_T,.[5i^l-^IV[SSl + ^^ 



TiiHH]- —Tr[SS] + —TI[T^'Tf,], (109) 



it is substituted by the form given in ( 102 ), see (see also ([7l|)). If compared with ( [7l| ) 
(the expression ( |102| ) is extended to include the meson fields S and T), ( |102| ) does not contain 
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the fields V, since in the CQM model the pions are not coupling directly to the meson fields 
and includes mass terms such as Ah,s,t- In the chiral Lagrangian approach for heavy meson 



states, see section |2.3.3| , the fundamental fields of the Lagrangian are the meson fields. CQM 
is a somehow more fundamental approach since it includes, together with meson fields, also 
the quark fields. When one adds in (71) the kinetic terms related to the S and T fields 
following a chiral Lagrangian approach, see [^], one must also subtract two mass shifting 
terms: 5msTr[SS] and (5m7^Tr[T'^T^], where 5ms = n^s ~ = ~ and 6mT = 
niT — rriH = At — Ah are defined in the mg oo limit. In the CQM model, that contains 
explicitly the heavy and light quark fields, 5ms and 5mT are substituted by As and A^^, i.e., 
the mass differences between the heavy meson masses Ms^t and the mass of the heavy quark 
involved mq; Ah comes in the H kinetic term. 

Ah is the free CQM parameter. We cannot deduce it from the model, but we can fix it 
by reasonable numerical values. On the other hand. As and can be computed once Ah is 
fixed. In the case of As, one only needs to solve (ll03| ). In the case of At, we will use some 
experimental information and a l/mq correction to the meson mass formula (we will discuss 
this point later on). 

The CQM expressions for IIh{Ah), Ils{As) and IlTiAT) are here given. They allow to 
calculate the renormalization constants Zh,s,t- 



UniAH) 

nT(AT) 



Zh' 



z 



-1 



h + {Ah + m)h{AH 
h + {As - m)l3{As) 
1 



{AH + m)^l4^ + I,{AH] 



(Ac 



m) 



BAh 
dh{As) 
OAs 



+ /3(Ac 



(110) 
(111) 

(112) 

(113) 
(114) 



and finally: 



3AI 



{Al 



+ (m + At) 



m 



(^ + At)^^^+/3(At) 



dh{AT) 
OAt 



dAi 



+ /i + 2Ar/3(Ar) 



+ /o + At/i 



(115) 



where = 1 GeV. The / integrals are given in the Appendix. 

At this point let us fix Ah and compute numerically As^t, the couplings G3, G4, and the 
renormalization constants Zh,s,t- 

We will consider everywhere in this work Hj-c^, Sren and Tren, but often, when notation is 
evident, we will drop the 'ren'. 

The Ah values will be taken in the range Ah = 0.3, 0.4, 0.5 GeV. A5 and G3 follow directly 
from eq. ( [1031) . From eqs. (Ill) and (|112| ) it is evident that As — Ah = if m ^ 0. Finally, 
At is obtained as follows. 

Take Mh and Mt to be the spin averaged masses related to the H and T multiplets 
respectively (a weighted average of the experimental masses of the particles in each doublet 
is taken. The weights are given by the number of polarization states that each particle can 
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assume according to its spin). We can write the following two 0{:p^) equations: 

Mh = mQ + AH + — (116) 
AL 

Mt = mQ + At + —■ (117) 

This couple of equations can be written both if mq = rric or mq = rrih. 

In the case of mq = rric we must use experimental information about the D2 and D\ states. 
As for D2, we have rn£)*^2AQo)o = 2458.9 ±2.0 MeV, r£).(2460)o = 23±5 MeV and ?7i£).(2460)± = 
2459 ±4 MeV, r£).(2460)± = 251^^7 MeV. These particles are identified with the 2+ states of the 
T multiplet |"^. As for 1)^(2420), experimentally it is found m£,*(242o)o = 2422.2 ± 1.8 MeV, 
r£).(242o)o = I8.9I35 MsV; this particle can be identified with the l"*" state of the T multiplet. 
We are ignoring a possible small mixing between the l"*" states belonging respectively to the 
S and T multiplets |35|. D2 and D\ states are quite broad states since their strong decays 



proceed in d wave, as pointed out in 2.3.3. 
From this analysis we obtain that: 

f A' — A' 1 

AT-Ai^ + -^^^ ^~ 470MeV. (118) 

rric 

If on the other hand we consider the mq = mfe case, experimental data on positive parity 
resonances show a bunch of states, not easily resolvable, having a mass Mb** = 5698 it 12 MeV 



and a width F = 128 it 18 MeV [49|. If we identify this mass with the T multiplet narrow states 
mass, we obtain: 

CA' - A' ^ 

At-Ah + — ^ 380MeV. (119) 

mj, 

Reasonable values of the heavy constituent masses are mj, = ms — 300 MeV and mc = m^ — 300 
MeV, where 300 MeV is the constituent mass discussed in 3.1.3| , while ms and mo are the 



experimental masses of the B and D mesons respectively. Solving simultaneously ( |118| ) and 
(|119|), one gets: 



At - Ai^ ~ 335 MeV. (120) 
The results for A5 and A^ as functions of Ah, are shown in Table ||; in Table || we list the 



Ah 


As 


At 


0.3 


0.545 


0.635 


0.4 


0.590 


0.735 


0.5 


0.641 


0.835 



Table 1: A values in (in GeV) 



CQM Gj and Zj values. 

Through A5, CQM predicts the following ms mass value (in literature these are the Dq, D^' 
states) : 

ms = 2165 ± 50MeV, (121) 
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A// 


I/G3 


Zh 






I/G4 


0.3 


0.16 


4.17 


1.84 


2.95 


0.15 


0.4 


0.22 


2.36 


1.14 


1.07 


0.26 


0.5 


0.345 


1.142 


0.63 


0.27 


0.66 



Table 2: Renormalization constants and couplings. A/f in GeV; G3, G4 in GeV ^, in 
GeV-i. 

where the central value is given in correspondence of A// = 0.4 GeV, while the upper and 
lower values are related to the remaining two A// values. This determination does not account 



for l/rric corrections (see eqs. ( |116[) and (117)). These states are very difficult to be observed 
experimentally since of their large width: from |^0[ one expects r(Do ~^ D^tt^) ~ 180 MeV 
and T{Dl' — > D*~^tt~) ~ 165 MeV. Theoretical predictions of ms available in literature are 



anyway larger, a typical value is ms = 2350 MeV. Very recent CLEO data |51] indicate, as 



the mass of the S multiplet, ms = 2461 MeV. This discrepancy of about 300 MeV between 
CQM and experimental data can be attributable to the absence of 0(l/mc) corrections in 
CQM calculations. Anytime we will need to use ms in applications, we will use both the CQM 
predicted value, for consistency, and the experimental one. 

3.3 extension to include p and ai resonances. 



The effective Lagrangian ( |93D for the CQM light sector can be extended to include p and 
ai resonances. The operative hypothesis needed is that of Vector-Meson-Dominance (VMD) 
(and of Axial- Meson-Dominance for oi). We will briefly sketch the VMD hypothesis: then the 
insertion of p in C^^ (analogously for ai) will turn out to be a simple step. 
Let us consider the electromagnetic form factor for vr"'': 

(7r+(p')|J^|vr+(p)) = {p + p')^,F^{t). (122) 

where t = {p' — p)'^. The best way to determine this form factor is to consider the process 



7P — > nvr"*" where 7 is a virtual photon coming from the scattering of an electron |52]. Now one 
does the hypothesis that F.^{t) is an analytic function in the variable t with a branching cut on 
the real positive axis. This hypothesis allows to write the following dispersion relation for 

^.(„.ir.,5«. (123) 

The lower integration bound is the threshold above which the form factor is different from zero. 
Observe that for t ^ 00 one has F7r(t) — > 0, since the probability of producing only two pions, 
when an enormous number of higher energy states are accessible, is extremely low. 

Let us suppose that i^7r(t) is dominated by the p^ resonance (VMD-hypothesis) in such a 
way that one can write, for the absorbitive part [5^: 



Im{F^{t')) oc S{t' - ml). (124) 

Then we have that: 



2 

m^p 



FAt) = (125) 



— t 



On the other hand, if we consider the diagram in fig. ^ |54|, where the j — p coupling is shown 
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Figure 6: 7 — p coupling diagram. 



we can write: 



O-ypfp 
m?p — 



(126) 



where the fp couphng with pions is an universal couphng, i.e., is the same with two nucleons, 
two p's, etc. fpS universahty is a consequence of electric charge conservation and of the 
complete p dominance. 



From eqs. (|125|J126D we conclude that the 7-/3 coupling constant is given by: 

,2 



9ip — ^ f ■ 
Jp 



(127) 



In a paper by N.M. Kroll, T.D. Lee and B. Zumino ||55[, it is shown that, at lower order in e, 

o 

TTl 

it makes sense to consider the interaction of an interpolating p field, -j^P^, with the photon: 

,2 

-p^'Ap, (128) 



Jp 



even if this term shows to be manifestly not gauge-invariant (one can prove that introducing 
an term in the interaction Lagrangian, this problem is solved). Equation (|128| ) can be 
considered as the coupling of the gauge field with the current: 

,2 



-Pp- 



(129) 



This identity between the p interpolating field and the Jp current makes sense only if is 
coupled to a conserved external current, as one can easily show writing the Lagrangian for the 
field as: 



£P = --{Fp,f--mypp+j-p, 



(130) 
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where an interaction j ■ p with an external current j (the p source) is included. If one takes the 
divergence of the equation of motion derivable from ( |130[ ), and if one considers that 5 • J = 0, 
then one has d ■ j = 0. The equation of motion is: 

{d'^ +m^)p^' = f. (131) 

Eq. (|131| ) produces the following relation between matrix elements: 



= (132) 



Therefore, at the level of matrix elements, we can see that the p and the photon sources coincide, 
modulo a factor j-. A problem arises if one observes that fp is not directly measurable at t = 0, 

due to the finite p mass: in p — > vrvr one measures /p for t = ^ 0. This means that one has 
to formulate the problematic hypothesis that fp is essentially constant in the (0, m^) interval 
(in this respect VMD is similar to PCAC). This problem is even stronger for ai, that has a 
mass larger than the p one by a factor of 3/2. 

Let us go back to CQM. We can couple the p interpolating field to the vector fermion light 
quark current: 

2 

TTl 

-j^Pt^X^X- (133) 
Jp 

The same thing can be made for ai, writing an analog interaction for the ai interpolating field. 
These interaction terms give Feynman rules for CQM vertices between light quark current -p, 
-ai: 

2 

TTl 

p vertex = i^Ye^ (134) 
Jp 

2 

Tn 

ai vertex = i^^f^^^e';, (135) 

where e and e' are the polarizations of p and ai respectively. 

The expression for the CQM effective Lagrangian must incorporate these results. Let 
us write using a notation mediated by the Hidden-Symmetry approach, which incorporates 



VMD [39|: 



+ x(«^^7m - T^)X 

+ x{A^lplb-ihpp^'-ip,-ihaa^'-i^-i^)x, (136) 



where: 



■J'{,x)pv — dpX,y djjX^ + [xp,Xi/], (137) 
describes the strength tensor for the p and oi fields. Moreover: 

Pp = ^"^^M 9V = 5.8 (138) 
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and, in an analogous way, we also write {nia — 1.26 GeV): 



_ .9A . 

On — I 



IJla 

9A = 9-5, 



(139) 



where p and a are Hermitian 3x3 matrices related to positive and negative parity light mesons. 
If we consider: 



gvfp 

9Afa ' 



(140) 

(141) 



then we are implementing VMD and AMD hypotheses and we are recovering (|134|135| ) vertices. 
Numerically one finds: 



hp ~ 



ha ~ 0.95. 



(142) 



For us fp = 0.152 GeV^, as emerges from and u> decays in e+e , and fa = 0.25 ± 
0.02 GeV^, as it comes out from r — > i/t-tttttt decay [^^. This result agrees with a determination 
made for fa by QCD sum rules g. Lattice QCD predicts fa = 0.30 ± 0.03 GeV^, [|9|. Since 
1/fa multiplies all amplitudes containing the ai meson, the uncertainty on fa will induce an 
uncertainty on normalizations for all the amplitudes involving the light axial meson. 



4 Strong Couplings 

4.1 Processes with one tt in the final state 

This is the first of two sections where the CQM model will be used to compute the coupling 
constants for the strong decays H — > Hit, S Hit, S Sn, T — > ffvr, T Sn and 
H — > H{p, ai), H ^ S{p, ai). A technique allowing to go beyond the soft pion limit hypothesis 
is also introduced. As showed previously, CQM incorporates a direct coupling of the pion to 
the light quark current. For strong processes with one pion in the final state, there is only one 
CQM diagram describing the decay of an heavy meson to another heavy meson and a pion. 
This diagram is represented in fig. ^. Different processes have different in and out heavy meson 
states and the soft pion limit hypothesis must be discussed case by case. The soft pion limit 
allows to simplify calculations, but it is a rough approximation if, e.g., transitions S — > Hit 
or T — > Htt are considered. In the exact chiral limit we can write, for the pion momentum, 
= {Qtt, 0, 0, q-jr) . In the heavy meson rest frame, v = (1, 0, 0, 0), we have v - q = Qt^. Moreover 
V ■ q = V ■ {p — p') = As — Ah ^ 0, where p and p' denote the momenta of the in and out 
mesons respectively and we have used the relation v ■ k = A, k being the residual momentum, 
(see 2.2), and A the mass difference between the heavy meson and the heavy quark there 



contained. Therefore, the soft pion limit is not very reliable for transitions such as 5 — > Htt, 
where — Ah ~ 140 — 190 MeV. We should observe that, if we adopt the soft pion limit 
hypothesis, the CQM diagram of fig. shows a very soft NG boson emitted from an internal 
line of a diagram, while we should have expeced an Adler zero of the emitting amplitude in 
this case. Anyway the CQM regularization scheme forces £^ in the loop to be quite close to m? 
and this saves the soft pion approximation (see discussion in |3^, pp. 175,176). 

In [^] the calculation of CQM amplitudes for the transitions H Htt and S Hit has 
been performed using in both cases the soft pion limit, which works as a good approximation 
only in the first case. In EO], a technique allowing to avoid the soft pion limit has been 
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Figure 7: CQM diagram for transitions meson — > meson + pion. 



introduced for the T — > Hit process. The same technique gives the possibiHty of improving 
also the S — > Htt calculation. Evidently, the soft pion limit is a good approximation for the 
S ^ Stt process Recent CLEO data [P] indicate that ms ^ uit ^ 2460 MeV so that the 
soft pion limit may be used also for the T — > Sir process [42]. In the next two subsections we 
will show how to compute the mentioned processes in the soft and not-soft pion hypothesis. 



4.1.1 H ^ Htt, the soft pion limit 

Let's consider the first term in the Lagrangian 

C = igTT[HH^-A^5], 

where the meson field H has been defined in (|6^). The transition V 
can therefore consider {D7r\iC\D*) and, using (|143D , we obtain: 



{DTT\i£\D*) = g 



Tr[P75 



(143) 

0~7r is allowed. We 
(144) 



where A has been expanded up to the first order in vr and the zeroth order in the expansion has 
been neglected. Observing that P = {H\P\\KC) and P* = ly KC\P*\H) can be made explicit 
using 



67|), the interaction term (144) reduces to: 



(145) 



where e represents the polarization of the 1~ state. As already observed, this interaction 
effective Lagrangian describes the coupling of the pion to the meson states, the fundamental 
fields at low energy, see fig. 

In the CQM model, where the fundamental fields are the heavy and light constituent quarks, 
the same coupling is modeled as a coupling of the pion to the light quark current. Figure |^ 
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Figure 8: In the low energy language of chiral Lagrangians, the pion is directly coupled to the meson fields. 



sliows a one loop CQM diagram containing two vertices meson-Q-x, see fig. ^ described by 
£hi introduced in (^) and one vertex pion-xx- This diagram can be computed as a standard 
loop diagram in the following way: 



TV[(7 • I + m)(-^7^75)(7 • ^ + m)75^7'^ea 



167r4 7 {i^ - rr{^Y{v i ^ ^n) 

Here is a legenda for the factors appearing in the preceding expression: 

• ( — 1) from the fermion loop 



(146) 



r" from the three quark propagators 

from the Feynman rules for the XX^ vertex, described in and for the two xHQ 
vertices, described in C^^; both carry a factor of (— i) 



Zh is due to the fact that the two meson fields coming in the loop integral must be the 
renormalized fields, ^fZu being the renormalization constant 



niH comes from the normalization conditions (6f:,67) 



Nc is the number of colours running in the loop {Nc = 3) 
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• {—j^jfij^) is the one pion expansion of the 7 • ^75 term contained in £'■'■. 



At this point one should calculate the trace in (|146| ) and the following integrals: 



Z„ 



z 









1 




167r4 j 




— m?)'^{v ■ 














167r4 j 




— m2)2(ti . 






( dH 








167r4 j 


(^2 


— m2)2(ti . 





(147) 
(148) 
(149) 



The technique for calculating Z with the proper time regularization procedure will be 
explained in detail in the next section, where the general case in which the two light propagators 
carry different momenta {i and i + q respectively) is analyzed. The soft pion approximation is 
achieved performing the ^ limit. The expression for Z is then equivalent to the expression 
for I4^(A}{) given in the Appendix. The Lorentz structures must be treated as in the following 
example: 



levr^ 



dH 



(150) 



In order to obtain Q and R, the contraction with v^v'^ and g^'^ in needed. In such a way one 
obtains integrals of the type described in the Appendix. 

A comparison between ( |145| ) and (|14(]| ) shows that, in order to obtain one has to extract 



from (146) the coefficient of e ■ q. The CQM expression for g is then the following: 



9 = Zh 



^/3(A//) -2(171 + ^/^h\ {h + AHhiAn)) - ^m^h{AH) 



where all the / integrals are listed in the Appendix. Numerically: 

g = 0.456 ± 0.040, 



(151) 



(152) 



where the central value corresponds to Ah = 0.4 GeV and the lower (higher) values, correspond 
to Ah = 0.3 GeV (A^ = 0.5 GeV). These values agree with what is found using the QCD 
sum rules method, g = 0.44 it 0.16 [^, A good aggreement is also found with relativistic 
quark models giving g ~ 0.40 [Sl| and g = 0.34 ||6^. CQM disagrees with the determination 
by Le Yaouanc and Becirevic |6S] according to which g = 1- 



The coupling constant g allows the determination of the hadronic width: 



6vr/2 



(153) 



Numerical values are described in Table ^. 

Table m, following |38|, includes also the B.R.'s predicted by CQM for the radiative decays. 
The CQM calculation of radiative processes proceeds without qualitatively new elements with 
respect to what explained so far. 

The soft pion approximation is also suitable for the 5 — > Stt process |E^. The meson 



interaction Lagrangian is given by the second term in eq. (72) 



C = ig'Tr[SS-fA-f5]. 



(154) 
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1 If^C S^'XT 






A rr — 5 HpV 


F.vn 






65.5 


70.1 


61.9 ±2.9 


D*^ 




34.5 


29.9 


38.1 ±2.9 


D*+ _ 




71.6 


71.7 


68.3 ± 1.4 


D*+ - 




28.0 


28.1 


30.6 ± 2.5 


D*+ _ 




0.4 


0.24 


-"-•-■--o.? 



Table 3: Theoretical and experimental B.R.'s for D*. The theoretical values are computed for 
Ah = 0.4,0.5 GeV. 

Once calculated the CQM loop diagram, a comparison -with the mesonic amplitude {STT\iC\S) 
gives: 



9' = Zs 



-m + -As 



{l2 + Ash{As))--m^hiAs) 



to -which correspond the numerical values: 



-0.13 ±0.04. 



(155) 



(156) 



Observe that g' can be obtained by g -with the simple exchange m — m. Ah As and 
Zh — > Zs- Indeed the loop integral for the process 1+ O+vr can be obtained by the loop 



integral for H Htt, shifting the 75 matrix contained in (146) on the right side of 70-. In so 
doing, one writes the expression that -would have -written for the S Sir process but having 



—7 • V, instead of 7 • in the projector. Expanding the trace in (146) in it's four terms different 
from zero, one realizes that this substitution in the projector is equivalent to an m ^ —m 
exchange. 

According to the CLEO data, mg ~ m^, therefore the soft pion limit is a good approxima- 
tion also for T Sir. The CQM approach gives [E5|: 



f = 2VZTZs[mV-T], 



(157) 



-where V and T are linear combinations of the Ij integrals discussed in the Appendix, while / 
is defined in (|73|). Numerically we find: 



/ = 0.91: 



hO.56 
-0.27- 



(158) 



4.1.2 T Htt, / 

Let's go back to the integral Z given in eq. ( |147| ) and evaluate it in the general case of pions 
bringing a momentum q.,^ ^ 0: 



Z 



iNr 



dH 



16^4 J (^2 _ j„2)((^ + g)2 _ j„2)(„ . ^ + ^ + ^g) ' 



(159) 



where the Feynman contour prescription is made explicit in the expression for the heavy quark 
propagator. Since the pion is the NG boson of the broken chiral symmetry, we put: 



= (g^,0,0,g^). 



(160) 



37 



Let's write ( |159D in the following way: 



Z = Hi+H2 



(£2 - m^){{i + g)2 - m2)(u . I + A - ie) 
- 2vri ndH-, or, (161) 

J (£2 - m.2)((^ + g)2 _ rr),2) ' V ^ 

where the proper time regularization is concerned and the Plemelij identity has been invoked 



^ P.V. — - — ±iTr6{x-xo), (162) 



X — [xq ± ie) X — xq 

in such a way that poles with negative real part lie in the upper complex £q plane. 
As a first step we compute Hi: 

d „ dx 

ill 



dm"^ J Jo {{i + qx)"^ - m'^){v ■ e + A) 

. d f „ dx 



d^' 



dm"^ J Jo {{£ + qx)l + m2)(i4 + A) 

d f ,4„ dx 



dm^j'^^'^Jo (4 + m2)(^4-iA(x))- '^^^^^ 

In the first equation the Feynman integral trick has been used while, in the second, the Eu- 
clidean rotation is performed. At this stage ^4 = —ign and we can shift £4 ^ £4 + ZQ'j^x as 
shown in the last equation in ( |163| ). Here A.{x) = A — q^^x. Now we must exponentiate 
the light quark propagator and cutoff the integrals according to the Shwinger's proper time 
prescription described in |3.1.3 We recall also the numerical factor and write {Nc = 3): 



167r^ dm2 Ji//^2 J J £4 — «A(x) 



167r3/2 J^ 



In the second equation we have used the following formula: 



+2 



/ ds-^ dxe"^ (^■)(l-erf(V^A(2;))). (164) 
Jl/AZ s^l^ Jo 



where: 



erf(z) = ^ / dte-^\ (166) 
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In a similar way let's work out H2: 



dm? J Jq {{i + qx)^ — m?) 



dm'^ J Jq ({i + qx)'^ + 



dm'^ J Jq (^I + m2) 



= 2vri-^ f dH Cdx /■'^'^^,e-(™^+^^-A^W) 
a"!'' J Jo Ji/A2 

where the same steps as for Hi have been followed. In the last equation is taken into 
account. Summing the expressions obtained for Hi and H2: 



16^3/2 



/ ds^j^ dxe'^ (^)(l + erf(V^A(a;))). (168) 
Ji/A2 s^/^ Jo 



The soft pion limit amounts to consider q.,^ — > 0, A(x) A and therefore Z /4(A), where 
/4(A), given in the Appendix, comes from the integration of (|147| ) . The CQM calculation of 
the process T — > Hit [Q, requires not only the calculation of Z, but also of Z^, Z^,^ and Z^^x 
since the CQM loop integral for T — > Hit is: 

i rrr^ u i^c ,4,Tr[(7 • ^ + m)7/,75(7 ■ + q) + m)^^{l + 7 • v)^,K] 



(169) 

Let's consider, for example, Z^y = CVf^v^, + Dq^j^q^ + E{Vfj_qu + Q/i^'i^) + Og^y. If we contract 
both members of this equation with the tensors v^j^Vy, q^qy,..we obtain a linear system of 
four simultaneous equations to be solved with respect to C, D, E, O, knowing the (integral) 
expressions v^v^Z^y, q^^q^Z^y, v^^q^Z^y and g^^Z^y. The matrix A, acting on the vector 
{C,D,E,0), contains powers of qj^ multiplied by numerical factors. The hadronic matrix 
element that one wants to compute, informs about the powers of q^ that should be eliminated 
in A. 

Let' recall now the eq. (170): 

C = ^TTlHTf^iiD^j ■ A)j5] + ^TilHT'^ii^ ■ DA^H] + h.c. (170) 



Using the same strategy followed in the calculation of g in 4J., we find that h' = hi + h2 is 
given by Q: 

h' = ^ZtZh{ — \i2 + AtZ{At) + ^ihiAn) - HAt)) 

+ P{Ri{AT),Si{AT),q^)}. (171) 

The polynomial P{Ri, Si, q-,^), given in the Appendix, is a sum of powers multiplied by some 
linear combinations of the integrals Jj and Z that we called Ri and Si. If we write down 
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the hadronic matrix elements for the processes T — > Htt (the following three), S Htt and 
H Hit (the last two), we understand that, e.g., in the CQM calculation of h' one should 
take account only of two powers of q^, therefore we neglect terms 0{q'^) {e.g., in A). 



{D+{p')7T-iq)\Dfip,7^)) 
{D*+{p',e)^~{q)\Df{p,r^)) 



{D*+{p',e)n-{q)\Dr{p,7^)) 

{D%p')7:+{q)\D+{p)) 
{D\p')7T+iq)\D*+ip,e)) 



rriT 



WsVue^m q +g 



{P ■ if 

rriT 



m 



W4- 



H 



ms 



(172) 
(173) 

(174) 

(175) 
(176) 



The processes involving the c quark are the most interesting since the charmed states, even 
the positive parity ones, are object of wide experimental investigation. The coupling constants 
gi, ..,g5 are connected to the coupling constants appearing at the level of meson Lagrangians 
by the following relations: 



91 
93 

94 

95 



52 = 2^mHmT 



h' 



2mHmT hi 
3 Kj; 

r 



(177) 
(178) 

(179) 
(180) 



where /' has been introduced in (73). 

Using recent data on the masses 
following widths: 



and the CQM value for h' = 0.65, one obtaines the 



r{Df 



r{Df 



*0 



^ D+TT-) 
D*+-K-) 

D*+TT-) 



4.59 X 10^^ MeV = 19.4 MeV 



h 



12 



1.33 X lO'^— MeV = 5.6 MeV 

h>2 

1.47 X 10^— MeV = 6.2 MeV. 

X 



h 



12 



(181) 
(182) 
(183) 



Due to the neutral pion decay channel, these widths should be multiplied by a factor of 1.5. 

According to the chiral Lagrangian for heavy mesons, the total width of the state is 
dominated by decays with only one pion in the final state. Therefore we can use the experimen- 
tal value for this width, 23 ± 5 MeV, to obtain h'^^^ = 0.51, in good agreement with h' = 0.65 



predicted by CQM. 

Equation (183) gives the total width for Df^ decaying to one pion, i.e., T = 9.3 MeV. This 
is only one half of the measured total width 18.913 5 MeV. This effect could be attributed to 
a mixing of r(l+) with 5'(1"'") or to strong 0{l/mc) corrections [65]. 
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4.1.3 S Htt, / 



Applying the technique developed in 4.1.2, we can compute the strong coupling /' describing 
S — > Htt. One finds: 



/' = V ZsZh 



Ri{As) - R2{As) - ^^^^ + m'Z{As) 



(184) 



In this computation only terms up to order have been considered. Numerically we find: 

/' = -0.76 ±0.13, (185) 

where the error is induced by the variation of Ah in the range 0.3, 0.4, 0.5 GeV. In this compu- 
tation the CQM value ms = 2.165ib0.05 GeV has been used (the error in this last determination 
doesn't affect much the /' numerical value). Anyway, recent CLEO data [^] give ms — 2.461 
GeV for the broad charmed state 5'(1'''). This discrepancy between the experimental value and 
the CQM determination, is most likely due to 0(l/mc) corrections. Since ms determines the 
value, the experimental value gives a different numerical result for /': 



/' = -0.56 ±0.11, 
consistent with a QCD sum rules determination |5C]: 

/' = -0.52 ±0.17. 



(186) 



(187) 



If one computes the process S Hit applying the soft pion limit hypothesis, the following 
numerical determination for /' is obtained |3^]: 



-0.85 ±0.02. 



(188) 



Neglecting mixing effects 5'(1"*") — r(l+) we can evaluate the width of D[ using the CLEO 
results for ms and the CQM value for /'. We get: 



T{D[ D*tt) = 240 MeV, 



(189) 



which probably accounts for the entire D[ width, due to the limited phase space. This result 
can be compared with the CLEO total width for the 1+ state, 290179^ ± 26 ± 36 MeV. 

4.2 Processes with p and ai in the final state 

In this section we will discuss the strong coupling constants HHp, HSp, HHai and HSai. 



These will turn out to be essential when the semileptonic decays B — > {p,ai)iv |39] are exam- 
ined. According to the notations introduced in |35], these couplings are parametrized in the 
following way: 



J^HSp = -iCTT{SHrPp) + ilJ^'^{SHa^^Hp)^^) 



(190) 
(191) 
(192) 
(193) 



The CQM approach for the determination of these constants amounts to a loop calculation 
of diagrams like the one in fig. ^, where the pion is substituted by a p or ai according to the 
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Feynman rules obtained in 3^, followed by a comparison with the hadronic matrix elements 
determined by the interaction Lagrangians above listed. 

Even if the strategy is the same as before, we have some new technical difficulties in dealing 
with p and ai, partly because of the polarizations of these states, partly because of the fact 
that p and oi are not massless. 

The first mentioned problem has influence on the expressions for the hadronic matrix ele- 
ments. Consider, for example, the transition H{1~) S{0~^)p. This process has contributions 
in s and d wave. The d wave contribution comes only from the Lagrangian term containing 
the factor p in (|191| ) while, the s wave contributions, come also from the Lagrangian term 
containing the ^ factor. With obvious notation, the matrix element is written in the following 
way: 

{p+{e,q)S{p')\H{7j,p)) = -ie*^^^(Sg^^+Dv^q^). (194) 

Applying the notations introduced in |3.3| , we can compute explicitly the traces in (|191D ob- 
taining: 

S = ^V^^T^(2C-4/i(i;.g)) (195) 
D = ^V^^I^(4/i). (196) 

Therefore to obtain S and D, a CQM calculation of C and p is required. Similar considerations 
apply for matrix elements containing the T state. 

Since p and ai are not massless particles, the loop integral Z, given in ( |197] ), must be 
examined in the case of a general q {q^ ^ 0) momentum. Once that Z is computed, the 
problem of determining Z^, Z^^,..is consequential. 

As we know, Z is given by: 

z=^r 

IGtt^ J {i^ -m'^){{£ + q)^ -m'^){vi + A + i€)' ^ ' 



Let's consider the following identity: 



1 1 



1 



(£2 _ m?) ((£ + q)2 _ ^2) q2 + 2£-q 
where q is, e.g., the p momentum. Then we can write: 



(^2 _ ^2) ((^ + qY _ ^2) 



(198) 



(199) 



Therefore, considering transitions (1) — > {2)p with (1) and (2) generic in and out meson states 
in the CQM diagram and calling x the mass value of p or ai, the following expression for Z 
easily follows: 



iNr f d^l 

z 



32x7r4 J (^2 _ ^2) 
where: 



{vi + Al)(^;' • ^ + f ) {v.t + A2){v'-i- f ) 



(200) 



/ v-q A1-A2 
V ■ V = = . (201) 
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Applying the notations given in the Appendix, Z can be also be written as: 

^ ^ /5(Ai,x/2,a;)- J5(A2,-a;/2,a;) 
2x 

where uj = v ■ v' and this equation is well defined for all a; values. At this point, we must care 
of integrals with several Lorentz indices: 

= ^lv^' + ^2v'^' (203) 
Z'"' = ^:ig^"' + VliV^'v'' + Vl^v'^'v"' + ^fi[vi'v"' + v'^'v''], (204) 

where the ilj are reported in the Appendix. The CQM results for A and (3 are: 

2 

TTl 

A = ^ Zui-^i + mZ) (205) 

2 

Tfl 

(3 = \/2 — ^Zj/[2mQi+mp02 + 2f^3-n4 + n5-m2Z]. (206) 
9VJp 

Here the functions Z, 0,j are computed with Ai = A2 = A//, x = nip, iv = mp/{2niB) where 
one takes the first l/mq correction to liJ = 0. Moreover: 

1^ 



rup 



V2gvL 



p 

2 



VZhZs i-ni - 2— + mZ] (207) 
V rUp J 



C = j^y^ZnZs {mp^2 + 21^3 + + n?Z) , (208) 

9vJp 

where Ai = A//, A2 = A5, x = nip and lo = (Ai — A2)/mp. 
The axial couplings of oi to -fT and S are here listed: 

Xa = -^^^/Z^s{-^i + 2n2— + mz\ (209) 



2 



= V2^^JZHZs{ma^2 + 21^3 - + ^^5 + ^^Z), (210) 



where now Ai = A//, A2 = A5, x = m-a and u = (Ai — A2) /ma- Moreover: 



= _^^ZH(m(z + 2^) (211) 



Ah, 


X = rUa, UJ = ma 


/{2mB). ] 


A 


= 0.60 GeV-i 


Xa 


/3 


= -0.86 


f^A 




= 0.16 GeV-^ 


fJ'A 


c 


= 0.01 


Ca 



Ca = ^^^ZH(-2mf]i + maf^2 + 21^3 + ^^4 + ^^5 + "^^^), (212) 

9 A] a 



: 0.85 X (0.25 GeVV/a) GeV"^ 
: -0.81 X (0.25 GeVV/a) 
: 0.23 X (0.25 GeVV/a) GeV"^ 
: 0.15 X (0.25 GeVV/a). 
As already mentioned, fa is the principal source of theoretical error on these results. An- 
other source of uncertainty is the variation of A// in the range A/f = 0.3 — 0.5 GeV (above we 
have used only the lS.11 = 0.4 GeV value). The latter produces a significative uncertainty only 
in the C,/?a,Ca determination since C = 0.01 ±0.19, /3a = -0.81^°;^^ and Ca = 0.15+°;J^ while, 
in the other cases, only a few percent variation is observed. A theoretical uncertainty of ±15% 
can be added to the constants A, /i, Aa, /^a- This follows, for example, from the calculation of 
A performed in |35] for the processes B* and D* — > D^. Other evaluations of the of A 

can be found in [S6| and p^]. 
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5 Semileptonic decays 



5.1 Semileptonic decays: leptonic constants 

CQM, as all constituent quark models, can only give model-dependent predictions that do 
not have the features of theoretical solidity achievable by the SVZ sum rules approach Q or 
by lattice QCD. Anyway, effective constituent quark models are sometimes easy to use tools 
allowing to evaluate the hadronic matrix elements. 

One of the main goals of the next future experimental physics is the high precision mea- 



surement of CKM matrix elements Vcb and V^;,; B physics experiments, such as Belle |68] and 
BaBar are already running in this direction. The semileptonic decays play a central role 
for the determination of CKM elements, therefore, a wide set of theoretical predictions for 
these processes, is a very actual need. 

In the next sections we will consider how CQM can help this need. The greatest part of 
the results we will obtain have already been determined through other approaches. 

The first target will be the determination of s and p wave semileptonic B decays to charmed 
mesons. The Isgur-Wise form factors governing the hadronic matrix elements relative to these 
processes will be calculated through CQM loop diagrams where an external weak current 
interacts with the heavy quark internal line, inducing a boost on the velocity of the heavy 
quark and possibly a change in its flavor. 

Next we will consider the semileptonic decays B — > p£u and B — > aiii'. An experimental 
determination by CLEO ||7^, ^ gives a branching ratio for the former amounting to: 

B{B° ^ p-£+u) = {2.5 ±OAt^Q ^± 0.5) x 10"^ (213) 

This process offers an important test for CQM. 

Experimental information are instead missing in the B — > aiiv channel. CQM has given 
the first prediction for the branching ratio of this process. A QCD sum rule prediction for the 
same process has been recently carried out in a paper by Aliev and Savci [0]. Discrepancies 
with this work will be discussed later on. 

A role of prominent importance in the determination of Vub is covered by the semileptonic 
channel B irtu and a precise measure of this process is one of the major goals of B- 
factories. CQM predicts a new contribution to the form factors governing B — > ii£v. This new 
contributions are due to CQM diagrams in which the vr is directly attached to the B loop. 
They influence the form factors inducing corrections between 10% and 30%, according to the 
moment carried away by the weak current. 

Let's start with the evaluation of the leptonic constants F and F~^ through diagrams of the 
kind of that in fig. ^. These constants will be useful to calculate the so called polar contribution 
to B ^ p, B ^ ai. Their definition is: 

{YAC\q^''^5Q\HiO-,v)) = i^T^v^F (214) 
(VAC|^7^Q|5(0+,w)) = i^/^v^'F^. (215) 

Computing the loop diagram given in fig. |9|, with the CQM rules one easily finds: 

F = ^ (216) 



= (217) 
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l+k 




Figure 9: CQM diagram for the lepton decay constants. 



where the renormahzation constants Z and the couphng constant G3 have been discussed in 
3.2. The calculation proceeds following the usual CQM strategy: from the loop integral one 



extracts the contribution and a comparison with the matrix elements 
F and . The numerical values are given in Table 0. 



and (|215[) gives 



Ah 


F 


F+ 


0.3 


0.33 


0.22 


0.4 


0.34 


0.24 


0.5 


0.37 


0.27 


values. 


Ah is 


expressed 



GeV3/2. 

Neglecting logarithmic corrections, F and F~^ are connected, in the mq 00 limit, to the 
leptonic decay constants fs and /+: 



{0\q7^75b\B{p)) 
{0\q^^b\Bo{p)) 

through the relations /b = F / ^Jms and 
one obtains: 



(218) 
(219) 



F^/V^B^- For example, if = 400 MeV, 



fs 



150 MeV 
100 MeV. 



(220) 
(221) 



The QCD sum rules analysis [^2| gives F = 0.30 it 0.05 GeV^/^, neglecting radiative cor- 
rections and even higher values, 0.4 — 0.5 GeV^/^, including corrections. Another QCD sum 
rules analysis |73] suggests F^ = 0.46 ±0.06 GeV^/^, which is somehow higher with respect to 
the CQM values. Lattice QCD WM gives fs = 170 ± 35 MeV. 
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5.2 b ^ c transitions 

We are interested in studying some weak heavy meson decays, i.e., decays with flavor changing 
of the heavy quark due to the presence of the charged weak current W. Let's focus on the 
b — > civ' transitions. Such processes are described by the following four fermion operator: 

O = ^^ct'^CI - 75)6i^7M(l - 75)1 (222) 

At the quark level one can calculate the matrix element: 

Ampq = {c£i>\0\b) = ^ Uc{pc)l^{l - lb)ub{pb)Uu{pu)-i ^{l - l->)ui{pi), (223) 

that is relevant only at very small distances and fails when, as is the case in the hadron world, 
quarks cannot be treated as asymptotic states. The hadronic amplitude Amp/j is: 

Amp^ = {D*lv\0\B) = ^^{D*\cr{l-l5)b\B){iu\u^^{l-^^)i\\AC). (224) 

The hadronic matrix term in Amp^ accounts for non perturbative QCD effects. It is 
possible to parameterize these effects through form factors depending on q^, where q is the 
momentum carried by the weak current. Strong interaction symmetries, and in particular 
HQET symmetries (in the limit mg oo), have influence on the properties of these form 
factors. Many times this influence acts in the sense of simplification, reducing the number of 
form factors needed to parameterize the matrix element. 



5.2.1 The Isgur-Wise function ^{uj) 

Let's consider an elastic scattering process B[v) B[v') mediated by the weak current V^. A 
model of what happens in the hadron is the following: due to the action of the weak current, 
the heavy quark is suddenly substituted by another heavy quark having the same flavor but 
different velocity v' . Light degrees of freedom must rearrange themselves in order to give rise 
to a -B meson moving with velocity v' . This rearranging process is mediated by an exchange 
of soft gluons (having momenta of order ~ Aqcd) with the heavy quark acting as the source 
of chromoelectric field. The larger is v' — v, the smaller is the probability to have an elastic 
transition: we have a suppression of the elastic form factor. In the oo limit, the form 

factor can only depend on v and v' in the scalar combination v ■ v' = uj. We can introduce a 
dimensionless probability function, £,{uj), that works as the form factor of the transition. This 
function is known, in HQET, as the Isgur-Wise function: 

{Biv^Kn^'h^lBiv)) = mBC{io){v + v'r, (225) 

where is due to a particular choice of the normalization of the external meson states. To 
convince oneself that there are no terms proportional to {v — v'), it is sufficient to multiply 
the matrix element in (|225| ) by {v — v')^ and remind from section |2.2| that 7 • vh^ = and 
hvl • V = h^. 

The interpretation of ^(w), as the probability for the elastic transition B(v) B{v'), 
suggests to assign a value of = 1 when lo = 1, i.e., when there is no change of velocity; 
smaller probability values (form factor suppression) are assigned when v ^ v' . 

We can examine this point in greater detail. The ^(o;) function describes the response of 
light degrees of freedom to the change of velocity of the static source of colour. This change 
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can also be accompanied by a flavor change of the heavy quark, li v = v' , the current causing 
the transition is therefore a symmetry current since light degrees of freedom do not resolve the 
flavor of the heavy quark (this is the flavor symmetry of HQET) . This symmetry current comes 
together with conserved charges, generators of the flavor symmetry, which are connected to the 
current by the well known relation: 



Q 



f'f 



J d^xht'{x)hy{x), (226) 



where f = f means that the current hasn't produced heavy quark flavor change. One can 
verify explicitly that the weak current responsible for the ui = 1 transition is a conserved 
current: 

df^Kv^K = 0, (227) 



where the property 7 • vh^ = hy and the equation of motion derived by (23) have been used. 



Meson states must be eigenstates of charge operators in such a way that: 

Qf'f\P{v)) = \P'iv)) (228) 
Qff\P{v)) = \P{v)), (229) 

where P denotes a pseudoscalar meson (let's simplify the discussion taking the case of J = 
mesons). If we write: 

{P'{v')\~K,rK\P{v)) = mp^{uj){v + v'r, (230) 

and if we consider v = v', take the /i = component, integrate with respect to x and use the 
(p8|), we get: 

C(l) = 1, (231) 
assuming the following meson state normalization: 

{M{p')\M{p)) = 2mMV^{2TTf6^{p - p'). (232) 

The normalization condition ^(w = 1) = 1 is particularly relevant not only because it 
is connected to the flavor symmetry of HQET, but also because it is not affected by l/mq 
corrections. The relevant corrections are ©(l/m-g). This is consequence of the Luke theorem 
[^], generalization of the Ademollo-Gatto theorem |7^], according to which, at a; = 1, there 
are no corrections to the hadronic matrix elements responsible for the semileptonic decays 
B D^v and B D*ii>. The leading corrections to the normalization of these matrix 
elements are, at the zero recoil point, of order l/rn^. Since Aq^j-j/m^ ~ 10%, zeroth order 
predictions in the l/rric expansion are, for u = 1, very accurate. Far from the zero recoil point 
= 1, corrections of order I/ttIq are suppressed by (w — 1) factors. 

The normalization condition ^(1) = 1 is therefore a good table test to understand if CQM 
allows for a correct calculation of the Isgur-Wise function. 

Let's go back to the weak transitions b ^ c. The decay 6 — > c is mediated by the left-handed 
current = C'j^{l — 75)6. This operator not only carries momentum, but also rotates the 
orientation of the the spin sq of the heavy quark during its decay. For an assigned value of total 
angular momentum of the light degrees of freedom si, the relative orientation of sq determines 
if the hadron in its final state is a D or D*. The heavy quark spin symmetry induces relations 
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Figure 10: Weak current insertion on the heavy quark propagator line. 



connecting the matrix elements for the decays B — > D£u, B D*tv in such a way that they 
could be expressed in terms of a unique universal form factor £,{00). 
Let's consider for example the following matrix element: 

{D{p')\V^^\B{p))=f+{q^){p + pr + f-{q^){p-pr. (233) 
where is the vectorial component of . In HQET one has : 

{D{p')\V^'\B{p)) = ^mBmDi{oo){v + v'Y- (234) 
Therefore we can see a relation between ( |233| ) and ( |234| ) : 

/± = Ccfe „ , =g(^), (235) 



where Cch is the correction that emerges from the matching with QCD, see 2.2.2 . In the 
following we will consider Cch = 1 for simplicity. In an analogous way one can consider the 
matrix elements describing the processes B D* , where and A^^ are both relevant. One 
finds: 

{D{v')\c-f^{l--f5)b\B{v)) = V^^^:^CH{v^ + v'^) (236) 
{D*{v',e)\c-f^{l--f5)b\B{v)) = ^mBmp, i{uj)[ie^,^pe*''v"'v^ 

- {l+u:)el + {e*-v)v'^]. (237) 

In order to determine ^(w) with CQM, one needs to calculate the loop diagram in fig. 
We will address the simple case in which Q = Q' = b. Let's observe that the momentum q is 
given by: 

q = {mbV + k) — {rrihv' + k'), (238) 
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and the momentum carried by the heavy quark having v' velocity is then: 

mftf + £ — g = ra\,v' + I + k' — 



(239) 



where we have used the relation (|238| ) . Using the CQM approach developed in so far, one can 
write the loop integral describing the diagram of fig. ^ in the following way: 



mHZH{-l)i'^i'' 



167r^ 



Tr 



75 



{{i - - 7v?){v' -{l + k' - k)){v ■ £) 



(240) 



where for simplicity we are considering as in and out states the pseudoscalar mesons of the H 
multiplet. We can perform the shift i — k ^ i and, observing that v ■ k = v' ■ k' = Ah, we can 
rewrite (24C) in the following way: 

Nr. 



3„-2 



mHZH{-l)i i 



Tr 



(£2 _ m'^)(v' .i^AH){v-i + Ah) 



(241) 



Once computed the trace, ^(w) is extracted by comparison with the transition amplitude (|225|) . 
The CQM expression for ^(w) is: 



^{oj) = Zh 



1+LJ 



hiAn) +{m + 



2A 



H 



I5{Ah,Ah,u;) 



(242) 



where the li integrals are listed in the Appendix. Let's observe that /5(A,A,1) = ^/3(A). 
Recalling the equation (p. 131), we then have 



e(i) 



1, 



(243) 



as expected. 

A very accurate determination of Vet can be obtained measuring the recoil spectrum of D* 
produced in semileptonic decays of -B. In particular, measuring the differential decay rate for 
B D*iv, one obtains an indirect measure of the product Vcb£,ioj). If we expand ^(w) around 
uj = 1 and if we suppose ^(1) = 1 (taking into account l/noQ corrections and bound state 
effects gives a value of ^(1) ~ 0.91 ||l^), then near u; = 1: 



(i-pV('^-i))K6, 



where p^jy^r, defined by: 



2 

Piw 



(1), 



(244) 



(245) 



works as a fit parameter allowing to extrapolate back to the Vet value that one could measure 
at zero recoil. 

The CQM numerical values for p'j^r are given in Table ^. What emerges is that CQM results 
are in a quite good agreement with those obtained with QCD sum rules: pj^r = 0.54— 1.0 ^f% . 
After an overall examination of the results obtained for the Isgur-Wise function in a series of 
quark models l7|, |7|, |l|, authors in @ find p'jy^, = 0.97 - 1.28. Lattice QCD indicates 
lower results, around = 0.64 p3[ |. 

In fig. 11 it is shown a plot of the (,{uj) given in ( |242| ) . 
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Figure 11: Isgur-Wise form factor vs. uj for different Ah values. 



5.2.2 Tii2{uj) and T'^I2{'^) form factors 

In this section we will consider the semileptonic decays of an H meson to S and T mesons. 
These decays are written as: 



B D**lu, 



(246) 



where D** can be either an S state, i.e., a charmed meson 0"*" or 1^ with si = 1/2, or a T 
state, i.e., a charmed meson 2+ or 1"^ with si = 3/2. There are two form factors describing 
respectively these decays: they are known as tx/2 and t^i2 |^^- The relevant matrix elements 
are: 



< Dl{v', e)|c7^(l - -iz)h\B{v) >= ^SmBmo* T3/2(a;) x 



[(u; + l)e>"-e, 



-a/3 



v'^v'^v'' 



(247) 



OJ] X 



ImBmn* 



(248) 
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{Do{v')\cj^{l-j5)b\B{v)) = ^mBTJiDo 2 t,/2{lo){v'^ - v^) (249) 
{Dl'{v', e)|c7^(l - 75)b\B{v)) = ri/2(a;){2 i e^^p^e*''v'^v'< 

+ 2[{l-u)el + {e* -vW^]}, (250) 

where I?| and are, respectively, the 2+ and l"*" states of the T multiplet, while Dq and D]' 
are the 0"'" and l"*" S states. 

In the preceding equations we ignore logarithmic corrections. The r form factors can be 
calculated by means of CQM loop diagrams, see again fig. 10. The strategy for writing down 
the loop integral is the same as before. We obtain the following results: 



Ti/2(^) = IH^s) - hi^H) + (Ah -As + m(l - u)) hiAn, A5, u;)] (251) 

' 2(1 — uj) 



and: 



T-3/2(^) 



^/Z]^Z^ r fhiAn) -hiAr)- {Ah- At) hi^H, At, uj) 



m 



L V 2{l-uj) 
UAh) + /3(At) + {Ah + Ar) hjAn, At,uj] 
2 (l + w) 

^ ^ 3S{Ah,At,u;)-{1-2u;) S{At,Ah,u;) 



2 {-I - U + LO^ + UJ^) 

+ {l-u;^)T{AH,AT,u;)-2{l-2ij)U{AH,AT,u;)\ , (252) 



where S, T, U are linear combinations of Ij integrals, see the Appendix. 

At a first glance it may seem that the r form factors are diverging for w — > 1, see for 
example Ti/2{^)- Let's take a; ~ 1 and neglect the term m(l — uj). Since in the heavy quark 
propagator in /3(A5) appears the velocity v' , one can write the term in square brackets of 
( |251| ) in the following way: 

fiH,., J-^fZ"'"!. . . (253) 



levr^y {i^ - m^){v ■ e + Ah){v' ■ e + As)' 

Define = — v'^ and, consequently, e ■ v = 1 — lo. The Ti/2{^^) limit for a; — > 1 is then: 

1 



lim 



1 Ti/2{uj) = lim Const. Ae ■ v, (254) 

1 ' e-*o e • V 



where A, in the e — > limit, is given by: 

~m^J - m^){v ■ i +''ah){v ■ £ + As) ^^""^^ ^^^^^ 
and amounts to (just contracting by v'^ and using the Appendix): 

A = ^^{Ah) + ^h{As) + ^(A^^ + As)h{AH, As, 1). (256) 

Since the phase space allowed for B decays to positive parity states is small, {uJmax = 1-33 
for Dl, D* and 

'-^max ~ 1.215 for , Dq), we can consider the following approximation: 

Tj{cv)^T,{l)x[l-p'^{u;-l)]. (257) 
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Piw 


n/2(l) 


Pi/2 


^3/2(1) 


^3/2 


0.3 


1 


0.72 


0.08 


0.8 


0.48 


1.4 


0.4 


1 


0.87 


0.09 


1.1 


0.56 


2.3 


0.5 


1 


1.14 


0.09 


2.7 


0.67 


3.0 



Table 5: Form factors and slopes. in GeV. 



Numerically we find the results given in Table [s] where we have a general scheme of all form 
factors calculated by CQM in so far. 

Let's compare CQM results with those obtained by other methods, see Table ^. As for T3/2, 
we have a good agreement with quark models. As for T1/2, we find a good agreement only with 



n/2(i) 


Pi/2 


^■3/2(1) 


^3/2 


Ref. 


0.09 


1.1 


0.56 


2.3 


CQM 


0.41 


1.0 


0.41 (mpiit) 


1.5 




8| 


0.25 


0.4 


0.28 


0.9 




7i 


0.31 


2.8 


0.31 


2.8 




81 


0.41 


1.4 


0.66 


1.9 




8§ 


0.059 


0.73 


0.515 


1.45 






0.225 


0.83 


0.54 


1.50 




5 B 



Table 6: Here Ah = 0.4 GeV. In the second paper in |73] a slightly higher determination for 
■^"1/2(1) is obtained using the SVZ method up to the next to leading order. 

In Table |^ we show the branching ratios for the s and d wave semileptonic decays of B 
to charmed mesons. Here Vcb = 0.038 tb = 1.62 psec. The results in Table |^ seem 
to contradict recent experimental claims that the broad resonances dominate and, therefore, 
they impose some further understanding [psl. A possible direction to look at, could be that 



of examining corrections, as pointed out in []89| . In the latter reference |8^], it is shown 
how in the relativistic quark model the corrections act in the sense of an approximately 
two-fold enhancement of the decay rates B DqIv and B Dl'ii' which, at the leading order 
of the heavy quark expansion, tend to be approximately one order of magnitude smaller than 
the decay rates B — > Dliv and B — > D2ii' (due to the Lorentz transformation properties of 
meson wave functions). 



Decay 
~B 
B 
B 
B 
B 
B 



Diu 
D*£u 



Ah-- 

3.0 

7.6 

0.03 

0.03 

0.27 

0.43 



0.3 



2.7 

6.9 

0.005 

0.008 

0.18 

0.34 



0.4 



Ah 
2.2 
5.9 
0.003 
0.0045 
0.13 
0.30 



0.5 



Exp. 

1.9 ±0.5 |9| 
4.68 ± 0.25 



0.74 ±0.16 [gO| 
< 0.85 



Table 7: Branching ratio (%) for the B semileptonic decays in charmed states via CQM. 
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5.2.3 The Bjorken sum rule 



A second important test table for CQM is the Bjorken sum rule. Let's introduce it briefly 
considering the weak decay of an heavy meson in which the quark b{v) is substituted by a 
quark c(v'). After the action of the weak current, light degrees of freedom must rearrange 
themselves to build the new charmed hadron state. There are many possible reconfigurations 
that can be assumed and we can associate an amplitude to each of them; the following sum 
rule holds: 

Am-p{b{v) c{v')) = Amp{B{v) Xc{v')). (258) 

Xc 

It means that we must sum over all the possible charmed final states. If there is a form factor 
suppression in the elastic channel, as we move far from co = 1, there is a compensating growth 
of the amplitudes involving excited states (as the positive parity ones). Indeed one can show 
that: 



2 

PlW 



i + EN/1(l)l' + M/W]- (259) 



This means that the compensation comes only from sg = 1"*^ states and p'j^r > \- Numer- 
ically we find that S and T states, (fc = 0), saturate almost completely the Bjorken sum rule 
for all Aff values. 



5.3 B piu, B aiiu 

The form factors parameterizing the semileptonic decays B piv and B aiiu are given by: 

<p+(e(A),p')l^7M(l-75)&|^«(p)> = ^^^^^W/3e>V^ 

TUB + rup 

- iel{mB + mp)Ai{q^) 

+ i{e*.q)^P±^A,{q^) (260) 
mB + rup 

2?77 



where: 



while for ai one has: 



<atieiX),p')\q%{l-j,)b\B\p)> = J^^^e^,^pe*yyP 

niB + ma 



- ie*p^{mB + ma)Vi{q^) 

+ ^(^*-«)|^^^^(«') (262) 
mB + ma 

+ i{e*-q)'^^q,[V^{q^)-V^{q^)l 
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where nia is the oi mass and: 



rriB + ma 



TUB - ma 



(263) 



Inia 2ma 

With this parameterization of the matrix elements [^], the following relations hold for q'^ = 0: 

^3(0) = Ao(0) (264) 
^3(0) = VoiO), (265) 

where q = {p — p')- 

5.3.1 Direct contributions 

For direct contributions we mean the contributions to form factors derived by CQM loop 
diagrams where the decaying meson couples directly to p or ai, see fig. f2l The Feynman rules 



needed for the computation of these diagrams are the usual ones, see also (134) and ( [1351) . 



The technique needed for computing the loop integrals has been developed in 4.2. The 



loop integral derived from fig. 12 



is: 



{-l){-i)y/ZHmH 



167r4 



reg Tr 



(7 • ^ + m)(-^7 • e)(7 ■ + q) + "i)(V, A)i^(-75) 



(£2 - m2)((£ + g)2 - m2){v i + Ah) 



(266) 



In this expression: 

• (—1) comes from the fermion loop 



(— i) comes from the vertex Q-Meson-x- The vertex with p (oi) doesn't introduce new 
z's since of (113^ 



-7 • e), where e is the p (ai) polarization, is the vertex described in (|134|). 



The CQM expressions for the form factors derived by the direct diagrams calculations, are the 
following (here we can find the Qi expressions introduced in | 

,2 



V^{q') 
Afiq') 



A^{q') 
A^{q') 



— — (^1 ~ i'^B + mp) 
Jp V ms 

Jp mB + mp 

[(m^ + mmpU))Z — tompQi — mpO,2 — 21^3- 

— — 2(Drig] 

[Z^ ( y o 9^6 V ^ ^ 

— \ mZ - ill - 2 [niB + mp) 

Jp \ mB \ mp) 



rrin 



TUB 



r^i I mpLU + 2m- 



m. 



mBj 



21^3 + 1^4 1 - 2 



m 



B 



+ 0.5 + 2^6 



ri 



mBmp 



7,2 ^1 , 
Z \ m — m h mnipUJ 

niB 



(267) 



(268) 
(269) 



(270) 
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Fi gur6 12: The CQM direct diagram. 



where: 



+ m„ — 



and: 



ri 



2mBmp 



2 2 2 

— q — rrip 



(271) 



(272) 



Z and are given in the Appendix. In the above expressions one must consider Ai = A//, 
A2 = Ai — mpu), X = nip. 

The calculation for the — > oi semileptonic transition proceeds in a similar way. We find: 



— i^i - mZ - — ) (mB + ma) 
Ja y ruB \ rua J 



(273) 



2m 



„-\J Zhviib ■ 

J a -mB + ma 

[(— + mmaU:)Z + 2mQi — (Dmafli + 
+{2mu) — ma)0,2 — 2Q3 — — — 2u}Qq] 



V,^{q') = !!i.[Il(mZ-n,-2^ + 2 — n,){mB + ma) 
Ja V mB V ma ma / 



(274) 
(275) 



Vo^iq') 



— -^JZumB r^i m^w + 2m — ^ ^ 2m + 

fa V \ m% mB J 

^2 (ma + 2m — ^— 2mu) \ + 29.^ + 9.4,(l- 2-^^— \ + $75+ 

V mBma J \ m^^J 



2Vte UJ 



mBma 



+ Z ( rr? + m — ^ mmaUJ 

* mB 
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^0 








^0 


F^(0) 


0.83 


0.69 


0.81 


0.33 


1.62 


1.13 


1.13 


1.13 


ap 


0.93 





0.87 


2.9 


1.13 


0.18 


1.3 


1.9 


bp 


0.02 





-0.17 


2.6 


0.12 


0.04 


3.8 


0.93 



Table 8: CQM direct diagrams contributions to the form factors governing the semileptonic 
decays B ^ p and i? — > ai. The -F^(O) values for the i? — > ai transition must be multiplied 
with the normalization 0.25 GeV^//a. The theoretical uncertainty is around ±15%. 



where now: 



2mBma 

111 

r[ = — -. (278) 



These results are amenable to a numerical analysis. Let's consider the following parametriza- 
tion: 




(279) 



for a generic direct form factor F^{q^) \ ap^bp have been obtained through a numerical analysis 
of the region going from to g'^ = 16 GeV^. Numerical values are listed in Table ^ 

The form factors describing the B ^ ai transition at = are proportional to the 
normalization factor (0.25 GeV^//a) (recall the problems in the determination of fa mentioned 
in |4.2| ). These parameters are also affected by a theoretical uncertainty of about 15%. 



5.3.2 Polar contributions 

The polar contributions to the form factors come from those CQM diagrams in which the weak 



current is coupled to B through an heavy meson intermediate state. In fig. 13 it is showed the 



related CQM loop diagram. The form factor will then have a typical polar behaviour: 

F^iq^) = -I^, (280) 

where mp is the mass of the intermediate virtual heavy meson state. This behaviour is certainly 
valid nearby the pole. Let's assume that it is valid all over the q^ range that we want to explore, 
i.e., also for small values. This hypothesis is a good one for the form factors Af,A2 (where 
the superscript P indicates that they are derived from the polar diagram), since they are 
numerically small, with respect to A^ , A^ , in the range of small values. Things are different 
for A^((^') and Vq {c^^, we shall come back on this point later. 



Using the strong couplings calculated in 4.2 and the leptonic decay constants obtained in 



5.1, we can calculate the different contributions to F^{0). As for the semileptonic transition 
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Figure 13: The CQM polar diagram. 



B ^ p, they are p9|]: 



Af{0) 



-V2gvXF 



m 



3/2 



B 



-(C — 2fiiOTnp 



ruBoimB + nip) 



Bo 



(281) 

(282) 
(283) 



where uj = mB/{2mp), while A,/i,/3, C and (7y have been defined and calculated in 4^. As for 
AQ{q'^), we have to impose the condition ( |264| ) ; a possible choice is: 



1 



(284) 



mpyJ2mB rn?B ~ ^ 

Let's discuss this equation in greater detail. The amplitude for the semileptonic process B ^ p 



can be written, see fig. 14, in the following way: 

Amp(S^/9) = (VACIA'^im^ 
F „ i 



m 



-{Hiq)p{k,e)\iC\H) 



B 



IvTLB q — 



Cie* ■ q), 



(285) 



where H is the state. We have used eq. ( pl4D and the interaction Lagrangian term multiplied 
by P in (|190| ). It is easy to find that: 



C = iV2gvP, 



(286) 
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H 



p(k,8) 



H(q) 



Figure 14: The polar diagram at the level of meson interactions. 



therefore: 



Amp{B ^p)= iV2gvP^= ^^-^ ^r^^ • (287) 



/niB (p' — m'^ 

Equation ( |287| ) should be compared with eq. (|26C 

i2mM3iQ')-A^iQ')V-^- (288) 

What then fohows is ( p8^ ). 

As for the semileptonic transition B ^ ai, one has: 

^^(0) = -V2gAXAF^^^^^^^ (289) 

1^1^(0) = ^^^^, (CA-2/.A^m,) (290) 
mB{mB + ma) 



VfiO) = -./2,W ^^"^f^"^"\ (291) 
where O = niB /{2ma). A similar reasoning as the one made before can be applied for Vq^((7^): 

Vi'iQ') = VfiO) + 9aPaF^ lo^y' „2 - (292) 

may IruB m^^ — 

One should note that, if we do the hypothesis of massless leptons in the final state, these 
and Aq form factors do not contribute to the process width and can therefore be ignored. 
Numerical results have been given in Table |^ together with the polar masses. A theoretical 
uncertainty of ±15% is estimated. In fig. |l^ we plot Ai and A2, while in fig. If: the form 
factors A, Vi and V2 are showed together. We don't plot V{q'^) in fig. 15 since the theoretical 
prediction for this value is affected by a large uncertainty. Plots do not account of the errors 
given in Tables. 
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^0 




vr 


^2 




mp 


-0.84 
5.3 


-0.11 
5.5 


-0.15 
5.5 


-0.019 


-1.48 
5.5 


-0.32 
5.3 


-0.57 
5.3 


0.07 



Table 9: Polar form factors for B ^ p and B ^ ai semileptonic decays. Polar masses have 
GeV dimensions. The -F^(O) values for the i? — > oi transition should be multiplied by the 
normalization factor (0.25 GeV'^/fa)- The theoretical uncertainty amounts to ±15%. 




Figure 15: Form factors Ai and A2 for the semileptonic decay B ~* p. 



5.3.3 Branching ratios and widths 

Using the numerical values given in the Tables, one can compute semileptonic branching ratios 
and widths. First of all let's compare CQM results for i? — > p with what is obtained by 
other methods. In Table |l^ we show the form factors obtained summing up direct and polar 
contributions: 

Here the CQM results are compared with those obtained by other approaches. In particular, 



I would comment on the method followed in |92]. The idea is that of describing the heavy 



meson field H as an effective interpolating field given by the operator product of the two 
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constituent quark projectors times a momentum space wave function, solution of the Bethe- 
Salpeter equation for the Richardson potential (having the confining behavior at large distances 
and the Coulombic one at short distances). With such a description of the heavy field, the 
calculation of the correlator (VAC|r{ Ji JJ2}|VAC), typical of potential models, where Ji and 
J2 are the effective interpolating currents of the light mesons in the final state. Mi and M2, while 
J is the current of the process, is substituted by the calculation of the amplitude (M1M2I J|//). 
This approach gives results for the semileptonic decays B — > pii' and B K*i^£~ that are 
in quite good agreement with those obtained by lattice and QCD sum rules (in the region 
where these latter methods are predictive). 



As one can see, the CQM result obtained for y ((7^) is affected by a large theoretical uncer- 
tainty since it turns out to be the sum of two numerical values almost identical in magnitude, 
but opposite in sign. Aside from this problem, the CQM values are generally in good agree- 
ment with those found by QCD sum rules; with respect to other approaches they show to be 
systematically higher. 
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CQM PMgl LCSR H SR g LLCSR || 

^'(O) -0.01 ±0.25 0.45 ± 0.11 0.34 ± 0.05 0.6 ± 0.2 0.35+^;^^ 

A^(0) 0.58 ±0.10 0.27 ± 0.06 0.26 ± 0.04 0.5 ± 0.1 0.27l[!:[J^ 

a({0) 0.66 ±0.12 0.26 ± 0.05 0.22 ± 0.03 0.4 ± 0.2 0.26t[!:[|| 

Ag(0) 0.33 ±0.05 0.29 ± 0.09 0.24 ± 0.02 0.30tg:[|^ 

Table 10: Form factors for the -B — > p transition at q'^ = 0. The CQM results are compared 
with those obtained by other theoretical approaches: a potential model (PM), light cone sum 
rules (LCSR), QCD sum rules (SR), lattice calculations with (SR), i.e. (LLCSR). The large 
theoretical uncertainty on y(0) is due to the strong cancellation between direct and polar 
contributions. 



Using Vub = 0.0032 and tb = 1.56 x 10"^^ s, we obtain for B piu: 

-- (2.5 ±0.8) X 10"^ 
: (4.4 ± 1.3) X lO'^ s"^ 



B{B° - 






To{B' - 


^P' 


-lu) 


r+{B' - 


-P' 


-lu) 


r_(5° - 


-p' 


-iu) 


(r+ + r_)(^0- 


^P' 





- (7.1 ±4.5) X lO'^ s"^ 
: (5.5 ±3.7) X lO'^ s"^ 

: (1.26 ± 0.38) X 10^ (294) 
where Tq, r+, r_ are referred to the three helicity states of p. This branching ratio is evidently 



consistent with the experimental one (213) 



The reported results may be computed for different V^h values just multiplying by |Kt;,/0.0032| 
and choosing a particular Vuh- The same holds for tb (sec) multiplying it by rB/(1.56 x 10~^^). 

A discussion on the theoretical uncertainty of these values is in order. The large error 
on 1/^(0) reflects in the determination of r+ and r_, but it doesn't affect Fq and has only a 
small effect on the branching ratio. Theoretical uncertainties for ^^'(0) and ^2('-*) ^'^^ probably 
related. The theoretical error on widths comes from the sum in quadrature of the error on 
y(0) and a ±15% error common to ^x(O) ^2('-')- 

Let's turn now to the semileptonic channel B a\iv. CQM predictions for it are the 
following: 

B(S° ^ aXiv) = (8.4 ± 1.6) X 10"^ 
To{B°^a+iu) = (4.0±0.7) X 10^ s"^ 
r+iB^^afiu) = (4.6±0.9) X 10^ s"^ 

r^{B^ afiu) = (0.98 ± 0.18) x 10^ s''^, (295) 

where Fq, F+, F_ label the three helicity states of ai. 

In the determination of these decay widths we have included only the uncertainty arising 
from fa- Lower values correspond to fa = 0.30 GeV^, while higher values correspond to 
/a = 0.25 GeV^. The theoretical errors arising from the form factors at = are difficult 
to estimate reliably and are not included in this analysis. We can guess that the theoretical 
uncertainty is larger at least by a factor of two. 

CQM predicts a branching ratio for the B^ a^iu decay higher than the branching ratio 
for B^ p^iu and this analysis shows that B^ — > aflu could be responsible for almost the 
50% of the semileptonic channel B X^ly] the B piv takes another 15%. 

The form factor A suffers for an analogous cancellation problem as V . The analysis of 
Aliev and Savci |pl|, based on the QCD sum rules method, gives a value for A that is five times 
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smaller than the CQM predicted one (with the sign reversed, but this is because of a different 
overall phase definition of the hadronic matrix elements). 



5.4 B 7t£u 



The semileptonic B ttIu channel is of crucial interest for the evaluation of the CKM matrix 
element Vub and it has been widely investigated in literature, see the reviews [^] and [p^]. 

CQM allows for a simple determination of the already known leading terms in the soft 
pion emission limit and of a new sub-leading contribution computed beyond the soft pion limit 
approximation, i.e., through the technique discussed in [4.1.2 . 

Let's then write the weak current matrix element for the semileptonic transition amplitude 



{^{q^)\V>^{q)\B{p)) 



9 " 



(296) 



where Fi(0) = Fo(0). 



5.4.1 The non derivative contribution 



The diagram in fig. 17 shows a non derivative coupling of the pion to the B meson and to the 



weak current. To obtain the contribution to the form factors deriving from the diagram in fig. 

up to the first order in vr, neglecting the term 



O, one has to expand x = ^Q^ defined in 3.1.2 



of order zero. In so far we have always considered x 
zero. 

Of course the diagram in fig 



q, truncating the expansion at order 



17 produces a result proportional to the leptonic decay 



constant of B, as one can easily see comparing fig. 17 with fig. |9[ 
On the other hand, the diagram in fig 



17 makes sense only for very small q-,^, i.e., in 



the limit in which almost the entire incoming meson momentum is brought away by the weak 
current. 



We therefore obtain 141 



ND 



fB 
fn 

/b 



(297) 
(298) 



where /b was defined in (|218|) . 



5.4.2 The polar contribution 



Consider now a diagram of the kind of that given in fig. 13 where, instead of p or ai , one must 
consider the pion attached, through a derivative coupling, to the light quarks. The intermediate 
meson state can be a 1~ state, belonging to the H multiplet, or a 0"^ state, belonging to the 
S multiplet. The latter case is less important since it produces only a small correction to the 
chiral limit. 

Consider the case in which the intermediate meson state is 1^. Having in mind fig. 
we expect a contribution proportional to gF, where g, the strong coupling constant HHtt, has 
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Figure 17: "Callan-Treiman" contribution to the form factors of the B tt semileptonic decay. 



33]: 



been calculated through CQM in in 4.1.1. The following contribution to Fi has been obtained 

Fg 



1 



If instead the polar meson is 0"*", one has: 



(299) 



— 



1 



(300) 



where h is the coupling constant HStt discussed in 4.1.2 and by B** we mean the 0^ state. 
The linear dependence of (pOOj) on makes this term not relevant in the chiral limit. 



The polar form factors ( |299 , |30C| ) are reliable near the poles, i.e., for q'^ ~ rns- We will 



discuss later about one way to extrapolate these results to a wider range. Numerically, if 
one uses the CQM results obtained in so far, one finds: 



0.52 ±0.01 
0.012 ±0.001. 



(301) 
(302) 



5.4.3 The direct contribution 

The contributions discussed above are well known in literature. CQM predicts the existence 
of a new contribution due to a diagram of the kind of that depicted in fig. ^ where the pion 
couples directly to B and to the weak current. This contribution differs from the non derivative 
one (ND), since it is derivative, and from the polar (which instead is derivative), since it doesn't 
contain couplings to resonances. 

The evaluation of these diagrams proceeds following the rules we have already discussed 
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many times. The result is: 



Zh_ 



[g^(C - mA) + mH B] 



(303) 



1 



m 



B 



j^iC-mA) 



m 



B 



where: 



A 
B 
C 



/tt V mu 
+ mu b(^ + 

-^(I3(Ah - q.) - HAh)) 
mA - m'^Z{AH) 

^{AHhiAn) - {Ah - qn)h{AH - Qn))- 



(304) 

(305) 
(306) 
(307) 



Here we have used the notation = ((77^,0,0,^7^), while Z is given in the Appendix. Observe 
that the dependence of ( |303| ) and ( |304| ) on q^ is due to the fact that: 



Numerically we find: 



qiT = ("T-l - q^)/2mB. 



Fl^'^iq' = 0) = F^''{q^ = 0) = 0.13 ± 0.05. 



(308) 



(309) 



The error made in the numeric evaluation is due to the variation of Ah in the range 0.3 — 0.5 
GeV. The direct diagram contribution induces a correction to the form factors that ranges from 
10% to 30% depending on q"^. 

Since the three contributions to Fq and Fi are independent, one can sum them up with the 
following result: 



Fj{q') 



fB 



(j + 1)/. 



1 — q'^ /m- 



(310) 



where mi = ruB*, mo = ms** and we have indicated the total form factor as F. F is the 
correct total form factor only when we are at the maximum value of q"^, i.e., q^ ~ = q^. 
The Fj form factors could be seen as the first terms in an expansion where higher terms 
are characterized by an increasing number of derivatives of the pion field. These terms are 
subleading only when q^ is high. 

Let's try to extrapolate to small values using the auxiliary functions Gj{q'^): 

F,{q')=F,iq')Gj{q'), 

where j £ (0, 1). This parameterization must satisfy the condition: 

GM) = 1, 

since ~ is the region where Fj are better approximating the Fj. 
Since Fi{0) = Fq{0), the following condition must hold: 



(311) 
(312) 



Fi(0)Gi(0) = Fo(0)Go(0). 



(313) 
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It is reasonable to assume that the corrections to Gj{q^) = 1 come from terms in which one 
has one more derivative of the pion field. Let's put: 



G,(,^) = l + -^ = l + i^, (314) 



where E^^ is the pion energy in the frame where B is at rest, while aj are free parameters. 
Moreover, since ( pl4D is equivalent to: 



the condition ( ^121) is automatically satisfied. Equation (|313| ) implies that ag and ai are related 
according to: 

(316) 



2Axai ruB \ Fo(0) 

niB V 2A;^aoy Fi(0)' 



It is not possible to fix the a constants only from experimental data, therefore we need to 
refer to other theoretical evaluations in literature. Quark models, @|-|@, predict Fo{0) = 
0.20 - 0.50, with the exception of |lOO|, which gives a value of Fo(0) = 0.09. QCD and fac- 



torization combined give Fq{0) = 0.33 |101], while chiral effective theory and HQET combined 



give Fo(0) = 0.38 ||. QCD sum rules give Fo(0) = 0.25 - 0.40 [|0|-[|0| and some lattice 



QCD computations give Fo(0) = 0.27 - 0.35 |l04|-pi)6i. 

Let's use, as an input parameter, the result obtained with QCD sum rules in [107], according 
to which Fq{0) = 0.30 ± 0.04. In so doing, one finds oq = 3.6. An oq > 1 suggests that the 
energy scale controlling the expansion (|314| ) is not A^ ~ 1 GeV, but larger. At this point one 



has to remind that the first term in (314) describes the situation in which q ~ Qq while, in the 



corrections depending on the pion momentum, we could have q.,^ values higher than expected, 
due to an effective expansion scale qA^ > A^. In such a way we can extend to small q^ values 
the range of validity of the form factor expressions. 

In Table |l^ we show, for some q"^ values, the two form factors Fi and Fq also containing 
the CQM correction. On the obtained results one can compute the error due to a variation 



of ao in the range (2.9,3.6,4.3). What can be read off from Table 11 is that CQM is in 
good agreement with the results obtained by other methods. This also means that the CQM 
deviations from the dominating contributions (the Callan-Treiman and the polar one), due to 
the direct contributions, are not strong enough to modify qualitatively the polar behaviour of 
the form factors predicted by the chiral effective theory. 
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14.9 GeV^ 


17.2 GeV^ 


20 GeV^ 


26.4 GeV^ 


CQM 










Ft 


1 c,8+U.28 

1 .Oo n p;o 

— U.oz 

— U.lo 


— U.oU 

0.62+0-Of 


2.96+!^-4^ 

— U.47 

0.65+0-°^ 

— U.IU 


13.78+n^^? 
0.83 ±0.01 


IS ro|] 










Ft 


0.83 
0.48 


0.96 
0.48 


1.19 
0.48 


3.14 
0.47 


GNS |98l 










Ft 
Ft 


0.82 
0.38 


1.05 
0.40 


1.45 
0.40 


2.31 
0.07 


LNS 


99|] 










pB-K 


0.53 
0.69 


0.57 
0.76 






Ball [ 


I03| 










pBn 


0.85 ± 0.15 
0.5 ±0.1 


1.1 ±0.2 
0.55 ±0.15 


1.6 
0.7 




Lattice (UKQCD) |106|| 










Ft 

pBn 


0.85 ± 0.20 
0.46 ±0.10 


1.10 ±0.27 
0.49 ±0.10 


1.72 ±0.50 
0.56 ±0.12 





Table 11: Form factors Fi and Fq predicted by CQM at high (near q^ax — 26.4GeV^) and 
comparison with other theoretical evaluations. The error quoted for CQM results is due to 
a variation of 20% in the parameter controlling the evolution from higher to lower values, 
where the results are less reliable. IS, GNS and LNS (from author names), are quark models, 
while the Ball's paper makes use of light cone QCD sum rules. 
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6 Appendix 



This Appendix is completely devoted to the Ij integrals used in the text and to their linear 
combinations arising in CQM applications. These integrals have been computed adopting the 



proper time regularization prescription discussed in 3.1.3. The analytical expressions here 



listed, have been numerically treated using Mathematica 3.0. Everywhere iVc = 3 is meant. 
/o(A) ' 



levr^ J (^v .k + A + ie 
Nc Z"^/^" ds 

16^3/2 



- ^Tii^H-'-A^vj ''''' 



^3(A) 



iA-, /■'■^s d^k 



levr^ 7 (A;2 _ • /c + A + ie) 

Ac /•^Z'^" _(is 

1/A2 S' 



4i e-^("^-^^) (1 + erf(A^)) (321) 



/4(A) 



16^3/2 

iNr d^k 



167r4 J [k"^ -m?Y(^vk + A + ie) 

t''^ e-(™'-^') [1 + erf (A^)] . (322) 
7l/A2 S^l^ 



167r3/2 



In these equations, 



t 



r(a,xo,xi)=/ (ite-*t"-\ (323) 



is the incomplete gamma function, while erf is the error function defined by: 



Let us now introduce: 



erf(z) = ^ / dxe-^\ (324) 



^.2 



^ A A ^ Ai (1-X) + A2X 

a(x, Ai, A2,uj) = — (325) 

^ ^ Vl + 2 (a;-l) x + 2 (1-u;) a;2 ^ ^ 
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/5(Ai,A2,Cj) 



iNr 



167r4 J {k"^ - m'^){v ■ k + Ai + i€){v' ■ k + A2 + ie) 



dx- 
6 



1 



X 



1 + 2x2(1 -a;) + 2x(w-l) 

ds a e-^(™'-"') (1 ^ erf(aV^)) + 

1/A2 



L 16^3/2 

6 



167r2 



ds e-^(-'-2-') "-^ 



1/A2 



(326) 



/6(Ai,A2,( 



167r4 



{v- k + Ai+ ie){v' ■ k + A2 + ie) 



lo + 2x2(1 - (j) + 2x(w - 1) 



dx- 



1 



167r3/2 7o 1 + 2x2(1 -w) + 2x(w-l) 



,(^2_,2) r + erf(^7V^)] • [1 + 2s{m^ - a^)] 
1/A2 s3/^ L 



(327) 



In the Ti/2,3/2 form factor determination, the following expressions are needed: 



S{Ai,A2,Uj) = Ai/3(A2)+L^ (/l + A2/3(A2))+Ai2/5(Ai,A2,Cj) 

T(Ai,A2,u;) = m2/5(Ai,A2,w) + /6(Ai,A2,cu) 

C/(Ai,A2,a;) = /i + A2 /3(A2) + Ai /3(Ai) + A2 Ai /5(Ai, A2, w). 



(328) 



Here are listed the integral Z, introduced in 4.1.1 , and all the auxiliary combinations used 
during the study of T — > Hn, S Htt processes: 



Z{A) 



167r4 



(A;2 — rin?)[{k + q)"^ — m?]{v ■ k + A + ie) 



f^/^" ds 



(329) 







where = ((^Tr, 0, 0, q-j^) is the pion four momentum and A(x) = A — xg^. Let us observe that, 
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in the soft pion limit q.,^ — > 0, one has Z[/S) /4(A). The auxihary combinations are: 



Ri{At) = m^Z{^T)-h{^H) (331) 
i?2(AT) = Ai.Z(AT) + (y + At) I2 (332) 

i?3(Ar) = ^{^tH^t) - I^hH^h)) (333) 

i?4(AT) = ^(/3(At)-/3(Ah)) (334) 

5i(At) = e - - A|/2 - Ai^Z(AT) (335) 

S2(At) = h + ^Th{^H)-m^h-m^/\TZ{^T) (336) 

2 

?77 

S3(At) = (/.(/i + A^f/3(Aj^)) + — (/3(A^^)-/3(At)) (337) 
54(At) = |^/i + ^(/3(Ah)-/3(At)) (338) 
S5(At) = ^^(A^f/3(AH)-/3(AT)), (339) 

where q-j^ = v ■ q. 

The following polynomial appears in the h' determination: 

PiRi, Si, q^) = --^ [SqlillmRi + 45i - 6^2) + 2ql{-nQmRi + 14Si + ^2 



+ 8S'3 + 48S'4) + 3q^{88mR3 + S3 - I6S4 - 32S5) + ISS's] , (340) 



while the following expressions emerge in the determination of / in |4.1.1 

V = i[(m2-A|)/4(AT)-/3(AT)-AT/2] (341) 

Y = ^ (342) 

r = ^[/i-y + (Af.-m2)/2 + AT(l- A^)/3(AT)-m2Ar/4(AT)]. (343) 

The processes mesons (p, oi) requires to compute Z in the case in which q^ = xv'^, co = v - v' , 
A2 = Ai — 2; w, X being the p{ai) mass: 



167r4 7 (A;2 _ f„2) p _^ ^^2 _ ^2] . + + ^g) 
h{Ai,x/2,oj) - h{A2,-x/2,uj) 



(344) 



2x 

The following fij expressions have been found in the determination of the strong couplings 



HHp,..., see 4.2, and have been used in the determination of the semileptonic form factors 
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B — > /o(ai), see 5.3.1 : 



K2 
K3 

^3 



4x 

/3(Ai)-3/3(A2) 



+ -[Ai/3(Ai)-A2/3(A2) 



xAi Ai[/3(Ai)- J3(A2)] ^ /3(x/2) - /3(-a;/2) 
-■^ H 1 



2 2x 4 

J3(-x/2) - J3(x/2) + L^[/3(Ai) - /3(A2)] [Ai - uJx/2]Z 

2a;(l-w2) 1-^2 

-/3(Ai) + /3(A2) - uj[h{-x/2) - h{x/2)] [x/2 - Aiuj]Z 



+ 



2x(l - ^2) 
2ujK4 -K2-K3 



2 ' 2(1 -a;2 



2(1-^2 

Acknowledgements 



2(1 -a;2) 


+ 


2(1-C^2)2 




-Ki 


+ 


3/^3 - Qu}Ki + K2{2uj'^ 


+ 1) 


2(1 -a;2) 


2(1-^2)2 






+ 


2i^4(2a;2 + 1) _ 3^(7^3 


+ i^3) 



2(l-a;2)2 



(345) 
(346) 

(347) 
(348) 
(349) 
(350) 
(351) 
(352) 
(353) 
(354) 
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